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Abstract: We study the non-perturbative effects of the global features of the 
configuration space for SU(2) gauge theory on the three-sphere. The strategy 
<^ is to reduce the full problem to an effective theory for the dynamics of the 

■ low-energy modes. By explicitly integrating out the high-energy modes, the 

i—i . one-loop correction to the effective hamiltonian is obtained. Imposing the 9 

! ! dependence through boundary conditions in configuration space incorporates 

£> \ the non-perturbative effects of the non-contractable loops in the full config- 

uration space. After this we obtain the glueball spectrum of the effective 
theory with a variational method. 

OO : 
O . 

^ ■ 1 Introduction 

^ ; One of the methods to gain understanding of non-perturbative phenomena in non-abelian 
gauge theory is by using the finite- volume approach Jj], ||. The mechanism of asymptotic 
freedom || renders the coupling constant small at high energies or, equivalently, at small 
distances. Taking the finite volume to be small results in a small coupling constant: one 

. £h ! can use perturbation theory. Gradually increasing the volume then allows one to study the 
onset of non-perturbative phenomena and to see, hopefully, the setting in of confinement. 
In this regime one can employ a hamiltonian formulation @ of the problem. The strategy 
in studying the dynamics of the Yang-Mills field is based on two observations. First, 
quantum field theory can in a sense be regarded as ordinary quantum mechanics, but with 
an infinite number of degrees of freedom. For the case of non-abelian gauge field theory, 
the configuration space of this quantum mechanical problem is complicated, which gives 
rise to non-perturbative behaviour. Second, coming from the perturbative regime, these 
complications first show up in the low-energy modes of the gauge field: one can capture 
the influence of the topology of the configuration space by studying the effective theory of 
the finite number of affected modes. 

Let us first reformulate quantum field theory as quantum mechanics. In the hamiltonian 
picture we are dealing with wave functionals in the configuration space. Speaking generally, 
each point in the configuration space corresponds to a configuration of the fields (p(x) . 
This means we have fields <p(x), conjugated momenta II(x), a hamiltonian H[(p, II] and 
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wave functional *&[(p]. If we decompose the fields (p(x) in orthogonal modes, we can 
reformulate the hamiltonian problem in terms of the generalized Fourier coefficients qk and 
their conjugate momenta Pk- This leads to a hamiltonian H[qk,Pk] and wave functionals 
^[(Z/c]. A typical hamiltonian would look like 



At small coupling, the modes will not interact and the potential will rise quadratically 
in all directions. As a consequence, the wave functions will be the familiar harmonic oscil- 
lator wave functions. One subsequently uses perturbation theory to take the interactions 
between this infinity of modes into account. What we have just described is the ordinary 
perturbative approach to quantum field theory. This perturbation scheme breaks down 
when at larger coupling the behaviour of the wave function for some of the modes is dic- 
tated by the properties of the configuration space, as explained below. We then have to 
replace these wave functions by functions that have the behaviour required by the config- 
uration space. In this regime the hamiltonian method is superior to the covariant path 
integral approach of Feynman. 

An essential feature of the non-perturbative behaviour is that the wave functional 
spreads out and becomes sensitive to the global features of the configuration space. This 
spreading out of the wave functional will occur first in those directions of the configuration 
space where the potential energy is lowest, i.e. in the direction of the low-energy modes of 
the gauge field (small ojk)- If the configuration space has non-contractable circles, the wave 
functionals are drastically affected by the topology when the support of the wave functional 
extends over the entire circle (i.e. bites in its own tail). This is what typically happens 
in non-abelian gauge theories. The Yang-Mills configuration space is the space of gauge 
orbits A/Q (A is the collection of gauge fields or connections, Q the group of local gauge 
transformations). We know from Singer || that the topology of this configuration space is 
highly non-trivial when Q is non-abelian. The configuration space also has a Riemannian 
geometry that can be made explicit, once explicit coordinates are chosen on A/Q. This 
geometry also influences the wave functionals, if the latter are no longer localized within 
regions much smaller than the inverse curvature of the field space. 

One of the most prominent properties of a non-abelian gauge theory, is the multiple 
vacuum structure. This means that the gauge degrees of freedom cannot be completely re- 
moved by imposing a simple condition like the Coulomb gauge. After gauge fixing there 
will remain Gribov copies, that is, gauge equivalent gauge field configurations that all sat- 
isfy the gauge fixing condition. Intimately related to this is the existence of instantons ||. 
They can be shown to describe tunnelling between gauge copies of the vacuum, that is, 
tunnelling from one Gribov copy to another. The configuration sitting at the top of the 
tunnelling path with lowest energy barrier is called the sphaleron |J: it is a saddle-point 
of the potential with one unstable mode. Let us suppose that, coming from the vacuum, 
the potential energy in the direction of the sphaleron(s) grows slowly. This means that at 
increasing volume the wave functionals around the different gauge copies of the vacuum 
will start to flow out over the instanton barriers and will develop a substantial overlap with 
each other. We then have the situation described above of probing the non-contractable 
loops in configuration space. 



H[qk,Pk] = ^2 \Pk + ^kQk) + interaction terms. 
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We impose gauge invariance by restricting the dynamics to a fundamental domain, that 
is, a set of gauge fields that is in one-to-one correspondence with the physical configuration 



space A/G |0]) 0, [[2). Explicitly, we first define A' to be the set of gauge field config- 
urations A with the property that || 9 -A|| > \\A\\ for all gauge transformations g. We will 
take this gauge field to represent the orbit on which it lies. One shows that A' is a convex 
subset of the set of transverse gauge fields, and that A' C f2, where Q is the Gribov re- 
gion which is characterized by the condition that the Faddeev-Popov operator is positive: 
FP(A) > 0. As A' is a convex subset of a linear vector space, it is topologically trivial. 
In order to get a one-to-one correspondence between A' and A/G, certain points on the 
boundary of A' (these are the points where the minimum of the norm of the gauge field 
on the orbit is degenerate) must be identified, after which we denote A' by A. Only after 
these boundary identifications, A has the non-trivial topological structure of the physical 
configuration space. If we restrict the dynamics to the fundamental domain A, the bound- 
ary identifications give rise to boundary conditions on the wave functional. See fig. |1] to get 
an idea of how the potential landscape might look in relation to the various spaces defined 
above. This figure displays a two-dimensional subspace (consisting of low-energy modes) 
of the configuration space. 

The use of the hamiltonian formulation is especially fruitful when the non-perturbative 
effects manifest themselves appreciably only in a small number of the modes of the gauge 
field. For this limited number of modes one defines an effective hamiltonian that takes 
the perturbative effect of all the other modes into account [[!]. By studying the resulting 
ordinary quantum mechanical problem, one can go beyond purely perturbative results . 

Apart from offering the possibility to go beyond perturbation theory, the finite- volume 
approach sketched above, which is sometimes called 'analytical', has the merit of being 
directly comparable to the numerical results of lattice calculations. For this the finite 
volume of the analytic approach must be chosen to be the finite volume of the lattice 
calculations. Putting a theory on a finite lattice necessarily means that it is put in a finite 
volume with certain boundary conditions. For the standard lattice geometry, this finite 
volume is cubic with periodic boundary conditions, i.e. a torus. 

To make the reduction of an infinite to a finite number of degrees of freedom, we use 
the Born-Oppenheimer approximation, which has its origin in the quantum mechanical 
treatment of molecules. In field theory one integrates out the fast modes from the path 
integral and obtains an effective hamiltonian in a finite number of slow modes. Using 
these methods in the intermediate-volume regime for SU(2) gauge theory on the torus 
resulted in good agreement with the lattice results || |I3| . In case of periodic boundary 
conditions on the torus and already at small volumes, one should take into account non- 
perturbative effects in the effective theory due to tunnelling through quantum induced 
barriers. These barriers are contributions to the potential coming from the integrating out 
of the fast modes. This tunnelling sets in much before the instanton effects play any role. 
Koller & Van Baal have shown that one can derive very accurate results, beyond the 
semiclassical analysis, by imposing non-trivial boundary conditions on the wave functional 
at the corresponding sphalerons. 

However one would like to extend these results to larger volumes in an attempt to get 
as close as possible to the confinement domain, and in this way get some intuition about 
how confinement sets in. The first thing to do, as argued above, is to include contributions 
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coming from instantons: one wants to take into account the sphaleron directions in config- 
uration space, because there the potential will show the strongest deviation from gaussian 
behaviour. This is not an easy task first of all because the instanton solutions on a torus 
are only known numerically [|l4j1 . 

To avoid this problem, we have shifted our attention to the case where the finite volume 
is S 3 , (the surface of) a three-sphere [15f|. Here the instantons are known analytically, and 
the tunnelling paths that connect gauge copies of the vacuum lie within the space of low- 
energy modes. Another simplification is that the structure of the perturbative vacuum is 
much simpler than in the case of the torus: there is no vacuum valley, that is, there is 
no continuous set of minima of the classical potential. A drawback of the three-sphere is 
that we lose the possibility to check against lattice results. Also, the approach to infinite- 
volume results will go as powers of 1/R as compared to the exponential behaviour in L 
for the torus fit) ). Thus, due to the intrinsic curvature of S 3 , it will be even less easy 
to derive infinite-volume results from our calculations. Nevertheless, within this model 
one can study non-perturbative phenomena like the effects of large gauge transformations 
and especially of the 9 angle on the glueball spectrum. A gauge transformation is called 
'large' if it cannot be continuously deformed into the identity. Typically, two Gribov copies 
of the vacuum are related by a large gauge transformation. The 9 angle, which will be 
properly defined in section |2|, is a free parameter of the theory that one has to allow when 
one implements gauge invariance under large gauge transformations. The conditions at 
the boundary of the fundamental domain will depend on the 9 angle, thus incorporating 
the non-perturbative effects of the non-contractable loops in the full configuration space. 
Previous research into gauge theory on S 3 [|17j did not include the study of the effects of 
large gauge transformations. 

This paper is set up as follows. In section ^| we define the effective model. In section |3| 



we explicitly integrate out the high-energy modes using a background field method [18 



This gives us the one-loop correction to the effective potential. In section ^ we construct 
a basis of functions that respect the (0-dependent) boundary conditions. Using this basis 
we perform a Rayleigh-Ritz analysis to approximate the spectra of both the lowest order 
hamiltonian and of the one-loop corrected hamiltonian. Summaries of these results ap- 
peared in respectively |19| and ||18|| . From the eigenvalues one can obtain the masses of the 
different excitations (the glueballs) in this model. Section |5] contains a discussion of the 
results. We argue that the effect of the instantons is large, but calculable. We determine 
the range of validity of the effective model and show that the results are in reasonable 
agreement with results on the torus. 



2 The effective model 
2.1 Introduction 

In this section we will define the effective model. First, the necessary machinery needed 
for the analysis on S 3 is developed. In particular, we will write down bases of functions for 
scalar and vector functions on S 3 . These bases of functions will be used in section |3] for the 
evaluation of various functional determinants. These functions will also allow us to isolate 
the space of low-energy modes on which the effective theory will be defined. After this we 
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will familiarize ourselves with some of the topological properties of SU(2) gauge theory on 
the three-sphere and we show how the instanton degrees of freedom are embedded in the 
space of slow modes. Then we will describe the derivation of the effective hamiltonian, and 
discuss the validity of the adiabatic approximation. The hamiltonian itself need to be sup- 
plemented with suitable boundary conditions in field space to obtain a well-posed quantum 
mechanical problem. These conditions emerge when imposing gauge invariance through 
restricting our problem to the fundamental domain. We will use dynamical considerations 
to argue that we have enough freedom to choose tractable boundary conditions. 



2.2 Low-energy modes 

We begin by introducing two framings on S 3 . These will lead to a basis of scalar functions. 
By taking tensor products of these scalar modes with eigenfunctions of the spin operator, 
we will construct vector functions on S 3 . By a similar method we will also construct 
su(2)-valued functions, where su(2) denotes the Lie algebra corresponding to the Lie group 
SU(2). After this we will write down the space of low-energy modes for SU(2) gauge theory 
on S 3 . 

Let r a be the Pauli matrices. We introduce the unit quaternions and their conjugates 

^ = °\ b y 

<r li = {%ir), a„ = (1,-if). (2) 
They satisfy the multiplication rules 

a i^v = V^a, °H°v = fjfr,<T a , (3) 



where we used the 't Hooft r\ symbols P0| , generalized slightly to include a component 
symmetric in /x and v for a = 0. The rf and ff form bases for respectively the self-dual 
and anti-self-dual four by four matrices. The rj symbols enjoy the following useful relations 



VpuVpu = VapV^p, Vu^Vup = Va^p, Vfr = Vlf3- (4) 

We can use rj and fj to define orthonormal framings of S 3 , which were motivated by the par- 
ticularly simple form of the instanton vector potentials in these framings (see section [2] 
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We embed S 3 in 1R 4 by considering the unit sphere parametrized by a unit vector n 
Using the scale-invariance of the classical hamiltonian, we can make the restriction to a 
sphere of radius R = 1. The R dependence can be reinstated on dimensional grounds. The 
framing for S 3 is obtained from the framing of H 4 by restricting in the following equation 
the four-index a to a three-index a (for a = one obtains the normal n M on S 
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The orthogonal matrix V that relates these two frames is given by 

Vj = ^ej M = ltT((n-a)a l (n-cj)a 3 ). (6) 
Note that e and e have opposite orientations. The parity operation V defined by 

V : (n ,n) -> (n , -n), (7) 
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interchanges the e and the e framing. 

A vector field can be written with respect to either framing (^). From now on we 
will, unless stated otherwise, use the framing and write 

K = A^. (8) 

Since we want vector fields to be tangent to S 3 we have that n^A^ = and the sum over i 
runs from 1 to 3. Latin indices refer to this framing, while Greek indices will refer to the 
embedding space 1R 4 . 

Each of the framings (||) defines a differential operator 



d i = e l d i = e l —— (9) 



'dx^ dx^ 

to which belong su(2) angular momentum operators, which for historical reasons will be 
denoted by L\ and L 2 : 



L\ = \d\ V 2 = \d\ (10) 

They are easily seen to commute and to satisfy the condition 

P = L\ = Ll (11) 

The manifold S 3 is isomorphic to the Lie group SU(2). The spatial symmetries of 
S 3 correspond to the left and right action of SU(2) on itself: so(4) = su(2) © su(2). In 
appendix [A] we show that L\ and L 2 correspond to the generators of the left and right 
symmetries on SU(2). 

A basis of scalar functions on S 3 is given by the set of eigenfunctions of the commuting 
operators {L 2 , Li z , L 2z }: 

\lm L m R ), i = 0, m L ,m R = (12) 

Since the laplacian on S 3 is given by <9,<9j = — 4L 2 , these modes are eigenfunctions of the 
spherical laplacian with eigenvalue — 41(1 + 1) and degeneracy (21 + l) 2 . See appendix |A] 
for other ways to arrive at this basis. 

To classify vector fields Ai on S 3 , we introduce a spin operator S by (S a A)i = —ie a ijA 
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This is an su(2) angular momentum operator that commutes with L\ and L 2 . We also 
define K = L + S , which is shorthand for = LSij + Sij. When there is no confusion 
possible, we will write L for L\. The fact that we single out L\ is related to the fact 
that we chose as the preferred framing for expressing vector fields. A basis of vector 
fields can be obtained by taking tensor products of the scalar functions ( |T2]) with the basis 
functions |lm s ) of the spin operator. In the standard way one obtains eigenfunctions of 
the set {K, L 2 }: 

\lm R ;km k ), 1 — 0,^,1,..., k — \l — 1\, . . . , I + 1. (13) 

The three modes with (I, k) = (0, 1) are the three vector fields e* . The three vector 
fields = V?e l correspond to the modes with (I, k) = (1, 0), as can be seen from V\Vl = 
itkimV^. To identify transversal and longitudinal modes, first note that d^A^ = diAi. Next 
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we introduce the operator Q via Qij = LiLj. Note that QA = for A in the Coulomb 
gauge. Using the identity 

(L-sf = L 2 - L- S -Q, (14) 

one concludes that the modes with k — I ± 1 are transverse, whereas the modes with k = I 
are purely longitudinal. 

To deal with su(2)-valued functions, we introduce the isospin operator T by T a = 
ad (|r ), where ad (X)(Y) = [X, Y\. This operator is yet another su(2) angular momentum 
operator that commutes with all operators defined before. By taking tensor products of 
the functions obtained so far with the basis functions |1 m t ) of T, we can construct scalar 
and vector modes that take their values in su(2). Introducing J = L + T we obtain for the 
su(2)-valued scalar functions (e.g. infinitesimal gauge transformations on S 3 ) 

\lm R ;jmj), l = 0,\,l,..., j = \l - 1|, . . . , I + 1. (15) 

For the vector functions (e.g. gauge fields on S d ) we define J = K + T and obtain 

\lm R ;k;jmj), 1 = 0, k = \l - 1|, ...,/ + 1, 

j = |A;-l|,...,A; + l. (16) 

— * — * 

Generalizing to the case where S and T correspond to respectively a spin-s and a spin-t 
representation, we obtain (suppressing the tjir dependence) 

\t; (ls)k; j rrij) , I = 0, |, 1, . . . , k = \l — s\, ...,/ + s, 

j = \k-t\,...,k + t. (17) 

The scalar modes are recovered by the choice s = 0, t = 1, the vector modes correspond 
to the choice s = 1, t = 1. 

An SU(2) gauge field on the three-sphere can be written as 

a, = a^ = A a 4y- ( 18 ) 

The field strength is given by Fij = diAj — djAi + [Ai, Aj] — 2eijkA k , where the last term 
is the so-called spin-connection: it is a consequence of the fact that we used a framing on 
the curved manifold S 3 . Reinstating the radius R would lead to a factor 1/R for this last 
term, which shows that this term vanishes in the limit of no curvature, i.e. for R — > oo. A 
gauge transformation g : S 3 x 1R — > SU(2) acts as follows: 

(*A) • = g- l A l9 + g-%g = g^D^g, (19) 

where D(A) denotes the covariant derivative in the fundamental representation. As usual, 
the field strength F transform under the adjoint representation: ( 9 F)ij = g~ 1 F i jg. In order 
to isolate the lowest energy levels, we examine the potential energy 

v ^ = ~^LM< F i)- ( 2 °) 
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Note the factor 2n 2 that we absorbed in the potential. We define the quadratic fluctuation 
operator Ai by 



V(A) = -— J ga triAMijAj) + O (A 3 ) . (21) 
This gives 

Mij = 2L 2 <% + 2K% - 4Q y . (22) 

Using S 2 = 2 and eq. ( |T4D one rewrites 

M l3 = (K 2 -P)\ y (23) 

Now focus on the modes ([16]). For I = 0, we have = 1 and Ai = A. The 9 eigenmodes are 
A^ 1 = cf with c" constant. For I = \, the modes with k = \ are pure gauge (.M = and 
Q / 0), whereas the modes with k — | are physical with .M = 9. For Z > 1, the modes 
with k = I - 1 have .M = (2/) 2 , the modes with k — I + 1 have .M = (2(Z + l)) 2 and the 
modes with k = I are again pure gauge. In particular the 9 modes with I = 1, k = have 
M. = A and are given by A" = c^W™ with c?^ constant. The lowest eigenspace of M. is 
thus 18 dimensional and given by 

A,(c, d) = (cf + rf^/) = (c^e; + ^) ^. (24) 

This is the space on which we will build the effective theory. 



2.3 Winding numbers and instantons 

Before writing down the instantons, we first introduce the notions of winding number and 
topological charge. The reader is referred to ]2T[ for more details. 

In the A = gauge, the set of gauge transformations Q consists of mappings g : S 3 — > 
SU(2) = S 3 . As a consequence, Q splits up in distinct classes that can be labelled by the 
winding number of the gauge transformations in that class: 

n[g\ = J sa %*tr ((g' 1 d i g)(g- 1 d J g)(g- 1 d k g)) . (25) 

Gauge transformations with zero winding number can be continuously deformed into the 
identity mapping, whereas mappings with non-zero winding number cannot. The latter 
mappings are sometimes referred to as 'large' or 'homotopically non-trivial' gauge trans- 
formations. 

Related to this, we define the standard Chern-Simons functional that measures the 
topological charge of a gauge field configuration: 

1 f 

Q[A] = J s3 e ijk tr (AiF jk - f^A^) . (26) 
The sign of n[g] was chosen such that we have Q[ 9 A] — Q[A] + n[g]. 
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In the hamiltonian picture, gauge invariance is implemented through Gauss' law, which 
implies invariance of the wave functional under small gauge transformations. For large 
gauge transformations, we have to allow for a possible 9 angle: 

iP[9A] = e in[9]e i)[A\. (27) 

The (anti-)instantons [0] in the framings (|5|), obtained from those on H 4 by interpreting 
the radius in 1R 4 as the exponential of the time t in the geometry S 3 x 1R, become (e and A 
are defined with respect to the framing for instantons and with respect to the framing 
for anti-instantons) 

Ao = , — r, A a = " A£ " , - (M + £ ;^, (28) 

2(1 + e-nY 2(1 + e-n) K J 

where 

2s 2 2sb u , 

U= l + V + s^ £ »= l + V + s^ S = Xe - (29) 

The instanton describes tunnelling from A = (Q=0) at t = — oo to A a = —o a (Q=l) 
at t = oo, over a potential barrier that is lowest when 6 M = 0. This configuration (with 
bp = 0, u — 1) corresponds to a sphaleron ||, i.e. the vector potential A a = := | a is a saddle 
point of the energy functional with one unstable mode, corresponding to the direction (u) 
of tunnelling. At t = oo, A a = — a a has zero energy and is a gauge copy of A a = by a 
gauge transformation g = n ■ a with winding number one, since 

n ■ ad a n ■ a = —a a . (30) 



We will be concentrating our attention to the c and d modes of eq. (|24|). These modes 
are degenerate in energy to lowest order with the modes that describe tunnelling through 
the sphaleron and " ant i- sphaleron" . The latter corresponds to the configuration with the 
minimal barrier height separating A = from its gauge copy by a gauge transformation 
g = n ■ a with winding number —1. The anti-sphaleron is actually a copy of the sphaleron 
under this gauge transformation, as can be seen from eq. (|28|), since 

n ■ oe a ^a a n ■ a = -e a ^a a . (31) 

The two-dimensional space containing the tunnelling paths through the sphalerons is con- 
sequently parametrized by u and v through 



A,(u,v) = (-uel-v<%)^=A i (u,vK, (32) 
Ai(u,v) = (-u8?-vV i a )^ = -u^ + vn-a^n-a. (33) 

The gauge transformation with winding number —1 is easily seen to map (u,v) = (w,0) 
into (u,v) = (0, 2 — w). In particular, as discussed above, it maps the sphaleron (1,0) to 
the anti-sphaleron (0, 1). Comparing with eq. ( p4[) shows that we obtain the (u,v) space 
from the (c, d) space by the choice c" = —u5f and d\ = —vSf. 
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2.4 Reduction to a quantum mechanical problem 

Let us start with the naive derivation of the hamiltonian for the (c, d) modes. From the 
lagrangian 

1 = -i? L = h I KK - if v(a) > (34) 

we obtain for the (c, d) space 

2tt 2 / • • n 2tt 2 

L = — 2 (c^ a + «) - -rV{c, d). (35) 

This leads to the hamiltonian 

f f d 2 d 2 \ 1 g 2 

K = lfe + « + f M) ' / = ^' (36) 



where the potential V(c, d) = V c j(c, d) can be obtained from eq. (|20|) : 

V(c,d) = \/(c) + V r (d) + |(tr(X)tr(r)-tr(Xr)), (37) 
V(c) = 2tr(X) + 6detc+ J (tr 2 (X) -tr(X 2 )) , (38) 



with the symmetric matrices X = cc T and K = dd T . 

The correct way to obtain the effective hamiltonian for the (c, d) modes must start from 
the full theory. In the case where the finite volume is the three-torus all zero momentum 
states become degenerate with the ground state at g — 0. This infinite degeneracy allows 
one to derive an effective hamiltonian for these states using degenerate hamiltonian per- 
turbation theory starting from the full hamiltonian 7i fll]. This full hamiltonian in the 
Coulomb gauge is given in Q. Proceeding to higher order in this way is complicated, 
mainly due to the non-abelian Coulomb Green function that occurs in the kinetic part of 



the hamiltonian, but it can be done 22 



For the case of the three-sphere however, the ground state at g = is non-degenerate. 
The (c, d) modes are in this respect not the analogue of the zero momentum modes on 
the torus, but are singled out just because of the fact that they are the slow modes of the 
system. To calculate the effective hamiltonian for these modes, we compute the effective 
action for the (c, d) degrees of freedom using a background field method. This explicit 
calculation will be performed in section [3]. 

We will use the hamiltonian formulation of the full problem for the discussion of the 
validity of restricting ourselves to the (c, d) space. Let x denote the 18 (c, d) modes and 
q all the transverse modes orthogonal to the (c, d) space. We use p x and p q to denote the 
conjugate momenta. The full hamiltonian in the Coulomb gauge [|J is a function of x, p x , 
q and p q . We define 

H\ x ] = Tt(x,p x = 0, q,p q ). (39) 
Consider the following decomposition of the full wave function 



* = £v (B) (*)xS(?), (40) 

n=l 
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where the functions |n) = (?) are chosen to be eigenstates of Tt[ x y. 

H [x] X ( 3(q) = V^(x) X $( q ). (41) 
The Schrodinger equation = E^J is equivalent to the following set of equations: 

(n\H(x,p s ,q,p q ) ( £ <p< m \x)\m)) = E^ n \x) Vn. (42) 



As an example we assume the following form for the full hamiltonian, 

n = -^ + jV dL {x)+n^{q,p q ). (43) 

This would be the form of the full hamiltonian if we neglect the complications of the 
non-abelian Coulomb Green function. The Schrodinger equation for this case leads to 

- { E VLv (m) (*) + (Vcl(^) + V^\x))^\x) = E^\x) Vn. (44) 
Here we introduced the covariant derivative V by 

V " m = lk 6nm + {n ^ m) = lk 6nm ~ A " m ' (45) 

The adiabatic approximation consists of truncating this equation to that is, we assume 
the transverse wave function to be in its ground state and we assume this ground state 
to decouple dynamically from the excited states. This approximation is valid if either the 
coefficients A l nm are small, or if we are in the region where the energy difference V <y2 \x) — 
VW(x) is large compared to the excitation energies of the states we are interested in. 

Returning to the general problem, the equations for the functions <p( n \x) will not just 
contain the functions A nm defined above, but more general matrix elements of the form 
{n\f(p x ,q,p g )\m). Still, if the excitation energies are small compared to V^(x) — V^>(x), 
the higher y9's will be small and we can restrict eq. ([42]) to ip^\ 



2.5 Boundary conditions 

We impose gauge invariance by restricting the dynamics to the fundamental domain A. 
The configuration space of our effective model consists of the intersection of A with the 
space of low-energy modes, the (c,d) space. This space was investigated in fll2|| , and the 
information obtained there on the boundary of the fundamental domain will be sufficient 
for our present purposes. To see this focus on the two-dimensional (u, v) plane (fig[j]), which 
contains the tunnelling paths. We have to provide boundary conditions at the boundary of 
the fundamental domain to obtain a well-defined quantum mechanical problem. At weak 
coupling, the potential energy at the boundary of the fundamental domain is higher than 
the energy E of the wave function: the wave function is localized around the perturbative 
vacuum c = d = and the boundary conditions are not felt. 

Increasing the coupling results in the spreading of the wave function over the config- 
uration space. We are interested in the regime where E is of the order of the sphaleron 
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energy: we will have a substantial flow of the wave function over the instanton barrier, but 
at the rest of the boundary the potential is still much higher than E. This means that 
at most parts of the boundary, the wave function will have decayed exponentially before 
reaching it. As a consequence, the boundary conditions imposed there will not have a large 
effect on the spectrum. By the same token, the precise location of the boundary in these 
regions is not important either. This gives us the freedom to choose tractable boundary 
conditions. At the sphalerons, however, the boundary conditions are fixed. Since the gauge 
transformation connecting the two sphalerons has winding number one, we have to set 

^(A(Sph, 0)) = e id V{A(Q, Sph)), (46) 

thus introducing the 9 angle. 

We define radial coordinates r c and r d by 

r e =[cfcf]*, r d (47) 

The sphaleron has radial coordinates (\/3, 0) and angular coordinates c" = — 5? (with 
cf = cf/r c ). It will be connected with the anti-sphaleron at (0, \/3). We will restrict the 
( r c r d) plane by r c < \/3, r d < a/3 and impose boundary conditions at the edges. When we 
come to the variational calculation, we will use basis functions of the form <f)(r c , r d )Y(c, d). 
As argued above, for values of the coupling constant at which our approximation will be 
valid, only the effect of the boundary conditions at the sphaleron will be felt. By imposing 
boundary conditions in the (r c , r d ) plane we pair up two submanifolds, of which only the 
sphaleron/anti-sphaleron need belong to the boundary of the fundamental domain. 
Consider the following decomposition of the full wave function 

1 oo 

* = ^E^M)xS](g), (48) 

' c' d n=l 



where q denotes all the modes orthogonal to the c and d modes. This is just eq. fl40|) with 
an extra factor r\r\ extracted for technical reasons (i.e. (p( n > «-> ■?/> (n ' ) r c ~ 4 r^ 4 ). Under the 
adiabatic approximation, we obtain a hamiltonian for ip = ift^ 1 ' given by 



d 2 . d 2 




H = -\ ^ + p- 12 b + 3l+^ + ^)+%4 (49) 



with Ac the laplacian in the angular coordinates. The extra term V^(c, d) that one might 
expect is contained within V(c,d). 

The boundary condition on ip follows directly from eq. (0), but care must be taken 
when imposing the condition on the normal derivative of ip. Matching along the sphaleron 
path across the boundary of the fundamental domain, we need to compensate for the 
curvature with the appropriate jacobian factor. We will make this more precise. 

Let us focus on the tunnelling path cf = — u 5f. Note that this path is equivalent to all 
paths c\ = —u S(a)f, with S an orthogonal matrix, due to the residual gauge symmetry. 
We will first remove this gauge symmetry to obtain a genuine one- dimensional tunnelling 
parameter. Introduce the following decomposition for c: 

c = S(a p )H(hi), S G SO(3), H T = H. (50) 
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If we write 



H = Y J h l H ll (51) 

8=0 

with tr(HiHj) = dij and with Hq oc i, then ho will play the role of the gauge invariant 
tunnelling parameter. The sphaleron is located at h = v^3, hi = 0, (i — 1, ... ,5). 

The isolation of the 18 (c, d) modes is appropriate close to the perturbative vacuum, 
because, as we have seen, they are the slow modes of the theory Close to the sphaleron how- 
ever, the only slow mode is the tunnelling mode w = h /\/3. To derive the correct boundary 
conditions, we should reduce the dynamics around the sphaleron to a one-dimensional tun- 
nelling problem, by integrating out all other modes. To this end, we consider the following 
decomposition of 

* = ^ n \w) X $ ] (q), (52) 

n 

where q not only denotes the non-(c, d) modes, but also the d modes and the hi modes 
with i = 1, . . . , 5. Note that \1/ does not depend on the gauge degrees of freedom a p . When 
imposing the boundary conditions, we want to relate the wave function at the sphaleron to 
the one at the anti-sphaleron. Consider the tranverse modes q within the (c, d) space at the 
sphaleron. The gauge transformation that maps the sphaleron to the anti-sphaleron will 
not map these modes on modes within the (c, d) space at the anti-sphaleron: transverse 
modes inside and outside the (c, d) space get mixed under the gauge transformation. The 
transverse wave function at the sphaleron restricted to the (c, d) modes does therefore not 
fit naturally to this transverse wave function at the anti-sphaleron. The full transverse 
wave functions do however map to each other. Another manifestation of this symmetry 
is the fact that the effective potential is symmetric around the sphaleron only when all 
transverse modes are integrated out. 

Before going to the effective hamiltonian for </?W [w) we must remove the residual gauge 
freedom. The gauge invariant wave function \1/ is independent of the coordinates a p . We 
will show that only after a suitable rescaling of the wave function, the laplacian takes its 
cartesian form with respect to the tunnelling parameter w. The boundary conditions at 
the sphaleron must be imposed on this rescaled wave function. 

Consider the laplacian for the c modes. If we write w M = (ho, . . . , h^, ax, 012, 013) the 
laplacian takes the form 



Here J = det 1 ^ 2 ^) and g is the metric which can be obtained from da 2 = tr(cc) = g^ v iXpU v . 
For \l/ independent of a, one can show that at the sphaleron path eq. (|53| ) reduces to 

= ^di (h%di$) . (54) 

To obtain an ordinary one-dimensional tunnelling problem, we must rescale the wave func- 
tion by w 3 / 2 oc ti^ 2 : 

3 

(f(w) = w~z(p(w). (55) 
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The tunnelling problem then looks as 
( f <9 2 



2dw 2 



+ V(w) 0(w) = Eip(w). (56) 



We can derive a similar equation for a function (p' of the parameter w' corresponding to 
tunnelling through the anti-sphaleron. When matching these tunnelling paths through 
the relation w' = 2 — w, boundary conditions must be imposed on (p and on its normal 
derivative at the sphaleron w = 1. Combining the factors w 3 ^ 2 and r -4 leads to the 
following boundary conditions on if): 

^(Sph,0) = e*V(0,Sph), (57) 
^M(Sph,0) = - e »^M(0,Sph). (58) 

3 The one-loop effective lagrangian 
3.1 Introduction 

In this section we will calculate the influence of the high-energy modes on the dynamics 



of the low-energy modes ||18|| . This is achieved by integrating out the high-energy modes 
in the path integral, which gives us the one-loop effective lagrangian for the low-energy 
modes. Denoting the low-energy modes by B(c,d), we will expand the one- loop effective 
lagrangian first up to the same order as the classical potential, i.e. up to fourth order in 
B and to second order in B. This allows us to do the renormalization of the lagrangian 
properly. 

Since we expect the physics to be sensitive to the details of the potential along the 
tunnelling path, we use an expansion of the effective potential around the sphaleron to 
construct a fit of the effective potential up to sixth order in the tunnelling parameter u. 
This allows us to write down the effective lagrangian including some fifth and sixth order 
terms in B so as to reproduce the behaviour along the tunneling path. 

It is our purpose to use the effective hamiltonian in a variational method to calculate 
the spectrum. Although it is possible to calculate the effective potential exactly along the 
tunnelling path, a polynomial approximation in B is much more useful, since this makes 
the analytic evaluation of the matrix elements feasible. 

The one-loop effective lagrangian will be given in terms of a path integral over the 
high-energy modes and over ghost fields. Using Feynman diagrams to expand these path 
integrals is the standard method to proceed, but due to the fact that the theory is defined 
on S 3 , the summation over the space components of a loop momentum looks rather different 
from the more familiar summation over the time component. To isolate the intricacies of 
the three-sphere, we will first assume B to be time independent. This allows us to perform 
the integrations over time components of loop momenta. Within this so-called adiabatic 
approximation the calculation of the effective potential as an expansion in B is then a 
purely algebraic problem. 

After this we return to the use of Feynman diagrams to evaluate the one-loop contri- 
bution to the kinetic term. We will neglect higher order time derivatives of B and also 
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terms of the form B n B 2 . We choose a renormalization scheme that makes the renormal- 
ized kinetic part assume the form of the classical kinetic term, replacing the bare by the 
renormalized coupling constant. The result of this section is then a finite, renormalized 
effective action. 



3.2 Gauge fixing 

We will impose the background gauge condition on the high-energy modes. Consider a 
general gauge field = (A , Ai) on S 3 , where A is the time component of the gauge field 
and Ai are the space components with respect to the framing e*. We will project out the 
background field B(c,d). Let Ps be the projector on the constant scalar modes, and let 
Py = P c + Pd be the projector on the (c, <i)-space. Explicitly we have: 



P#l> = P ^=wl^ (59) 
(PcA), = (^oH = i/ S 3^ ( 6 °) 
{P*A)< = {P^)=(^-jA k vi]vi (61) 



,27r 2 7s 3 

We define the background field B and the quantum field Q by respectively 

5 M = (PA)„ = (P S A , (P v A) t ) , Q„ = A li - Bp. (62) 

We define the gauge fixing function % = (1 — P s ) D^(PA)A^ + PsA . We use x to impose 
the background gauge condition: 

X = 0^l B \ s = ° (63) 
Introducing the Faddeev-Popov determinant 

A[A] = (jDgS\x( g A)]y\ (64) 
and performing the standard manipulations with the partition function leads to 



J DB k D'Q^[A] exp i j tr (±F%(B + Q) + ~(D„(B)Qrf 



(65) 



The primed integration means that we have excluded the / = modes from the integration 
over the scalar field Qo and the 1 = 0, k = 1 and 1 = 1, k = modes (the (c, d) modes) 
from the integration over the vector field Qk- 

We now focus on the Faddeev-Popov determinant (|64[) . Under an infinitesimal gauge 
transformation A the change in A^ is given by D fJi (A)A and the change in the gauge fixing 
function by 

= (1 - P s ) {D,(PA)D,(A)A + [(PD(A)A),, A,}} + P S D (A)A 
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= (1 - P 8 ) {D,{B)D,{A)A - [B„ (PD(A)AU 

-[Q lt ,(PD(A)A) ll ]} + P s D (A)A 
= (1 - P s ) {D^B){1 - P)D,(A)A + d fl (PD(A)A) fl 

-[Q„(PD(A)A) fl ]} + P S D (A)A 
= (1 - P 5 ) {^(5)(1 - P)D M (A)A - [Q„ (PD(A)A),)} 

+P S D (A)A. (66) 

The Faddeev- Popov determinant is given by A [A] = det(Sx/SA). Introducing ghost fields 
i/j and ip, this determinant can be written as a fermionic path integral. We will evaluate 
this path integral in the quadratic approximation, that is, we throw away terms of order 
three and higher in the quantum fields ip and Q^. After this the integration over the I = 
modes of if) can be performed to give an irrelevant constant and we obtain 



A [A] oc J Dty Dty exp (3 J tr (^(^)(1 - P)D lx (B)ij; 



(67) 



Substituting this in the partition function with j3 = —l/gfi, expanding the classical eu- 
clidean action up to second order in Q, and choosing the Feynman gauge £ = 1 we obtain 



DB k D'Q^ D'i) D'lp exp 



1/tr U{-D,{B){l-P)D,{B)}^ 
9o J 



- ±F%{B) - 2{D ll F ilv ){B)Q v + Q„W I ^(B)Q I 



with 



W 00 



-D 2 P {B) 

-Wio = -2 ad (Bi) 

-2 ad (F ij (B))-(D 2 p (B)) ij + 25 ij 



(6J 



(69) 



Remember that the covariant derivative Di(B) acting on vectors (or tensors) gives extra 
connection terms (due to S 3 being a curved manifold), e.g. 



diCj + [Bi, Cj] — SijkCk 
(-2iLi + iB?T a 



(70) 
(71) 

The action contains a term J U Q U with J u = (D fl F /JiV )(B). Since B need not satisfy 
the equations of motion, this term does not vanish. Upon expanding the path integral in 
Feynman diagrams, this term will give rise to extra diagrams, where J acts as a source. 
In appendix [FJ we will show that the presence of J will only contribute to terms in the 
effective lagrangian that we will consider to give only small corrections: they are at least 
of order c 2 (P or c A d 2 . 

Dropping for the time being the term linear in Q p , we are left with a gaussian integration 
over the fields tjj and from which we extract the effective action 



sf[B] 



[ dr()C + V) 
Jo 

S$[B] - lndet {-D„{B){1 - P)D fl (B)) + \ lndet (W t 



(72) 

Although we will not use the primed notation to denote it, the determinant does exclude 
the same modes that were not integrated over in the path integral. 
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3.3 The effective potential 



For computing the effective potential, B is considered to be independent of time. The 
assumption B = directly implies Woi = 0. This results in a factorization of the integral 
over Q M , and we obtain the one-loop contribution to the effective action 

S$P[B,B = 0] = -lndet(-D M (fl)(l-P)£V(fl)) 

+ | lndet (—D fl (B)D fl (B)) + §lndet (Wij) . (73) 

We have to calculate these functional determinants for a general background B(c, d). 



3.3.1 The operators 

We start by expressing the various operators of eq. flT3|) in terms of c and d. For the scalar 
operator Woo = —D 2 (B) we find 

- D^B)D^(B) = -d 2 - D t (B)D t (B). (74) 
As it acts on scalar functions, we have Di(B) = —2iLi + iBfT" 1 and we obtain 

- Di(B)Di(B) = AL 2 — AB?LiT a + B1B b T a T b 

= AL 2 - Ac1L u T a - Ad^L 2i T a 

+ (cfc{ + d«d\ + c«d b m Vr + d a m Vrci) T a T b . (75) 

The other scalar operator (the ghost operator) can also be written as 

- D,{B){1 - P)D,{B) = -d 2 - A(£)(l - Pv)D i (B), (76) 

where use was made of the fact that the operator does not act on constant functions. 
Carefully analysing the extra term allows us to write 

A(B)iVA(B) = P^D i (B){p v ) ij D j {B)p 1 L 

= -\{ (c?c b - ctc)L l3 L u ) + (d*d\ - d?d b L 2j L 2l ) } T a T h Pt (77) 
The vector operator W^ is more complicated due to the extra connection terms. We write 

Wij = -d 2 5 tJ + (78) 
Using the spin operator S we can suppress the spin-indices and write 

W{B) = 2L 2 + 2K 2 - Ac a k L lk T a - Ad a k L 2k T a + (-6cg + 2d a m V k m ) S k T a 

+ ( ^ijk^abcCiCj 2SijkB a b c C i d m Vj + Elmn^abcdi d m V k ^ SfcT 

+ {c a A + d a k d\ + c a k d b m V k m + d a m V k m c b k ) T a T b . (79) 
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3.3.2 Reduction to spatial parts 

As we have seen, the operators we are interested in are of the form 

A{B) = -dl-dl + A{B), (80) 
where we introduced a laplacian term for an e-dimensional torus of size L that we attached 



to our space S 3 [23|. This allows us to neatly perform the dimensional regularization. The 
scale L should of course be chosen proportional to the radius R of the three-sphere. A 
precise choice of L would fix our regularization procedure completely, and different choices 
will be related through finite renormalizations fixing the relations between the associated 
A-parameters. 

Suppose that the spectrum of A is {Aj}. The spectrum of —d\ is {k^} with k e = 2nn/L. 
If we also take the time periodic with period T, the spectrum of — <9q is {k^} with k = 
2vm/T. Since [d ,A(B)] = [d B ,A(B)] = 0, the spectrum {A;} of A{B) follows trivially 

Consider the exponential integral 

r°° ds 00 (— Tl n r n 

E l (x)= ^e- s = - 7 -ln(*)-£^L^. (81) 



This implies 

ds 



x s ~! nn\ 



-As 



-ln(A) -MS) -7 + 0(A<J). (82) 

<5 S 

Taking the limit 5 j in this integral, we can write, up to an irrelevant constant, 
lndet(A) = ^ln(Ai) 



EEE 



°° ^ e - S (fc 2 +fc|+A s ) 
S 
£ roc 



= r ,fa '" , " , * r ( e '^) <83) 

where we replaced the summations over k and k e by integrals. The operator A can be 
expressed in terms of the angular momentum operators defined above, the functions V/ 
and the constants c" and df. To take the trace, we need a basis of functions (see section 0). 
For the scalar operators (the ghost operator and W 00 ), we can use \lm L m R ) ® |lm t ) or 
equivalently \l iur; j rrij) , where mi, rriR and rrij correspond to the z-components of Li, L 2 
and J = Li + T respectively. For the vector operator W{j, we can use \l tur] k nik) ® \ 1 m t ) 
or \l rriR] k; j rrij) . Note that the c and d modes correspond to the vector modes with 
(l,k) = (0,1) and (l,k) = (1,0) respectively. For the scalar operators, the trace must 
not be taken over the I = functions, whereas for the vector operator the trace must not 
include the c and d modes. For the case of the ghost operator, the operator (1 — P) (see 
eq. (|76|) ) causes some intermediate vector modes to be projected to zero. 
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3.3.3 Exact determinants 

Using the appropriate basis, we can calculate the determinants exactly for the vacuum 
u = (5 = 0), and for the sphaleron configuration u = 1 (c" = —Sf and d = 0). After 
integrating over s, the final summation over / is expressed in terms of the function ((s, a), 
which is defined by 

OO J 

tf a ' a ) = £(5fcM^w Re (*)>i ( 84 ) 



fc=2 



After analytic continuation we have 



OO 



C(a,o)=X; i %(-or(U2a + 2m)-l) ) S ^|,-± (85) 



m=0 ' Ul 



where Cr denotes the Riemann ^-function and (s) m is Pochhammer's symbol. If s ap- 
proaches one of the poles, this expansion can be used to split off the divergent term: the 
remainder of the series is denoted with (j?(s,a). As an example we treat W(B = 0) = 
-8q + 2L 2 + IK 2 . From eq. @ we obtain 

lndet(Wy = -^=(^=) £ r(-i-|)tr((2L 2 + 2i? 2 )^f) 





' L 
















f L 




K 2^ 



) r(-|-i){3C(-i-| -i) + 3C(-| -l-i) 

+6 (Cji(-3 - e) - 9) - 6 (C«(-l - e) - 3)} . (86) 
The summation over (/, k) is as follows: 

OO 

^ = ^ (Wta + ^ra+l^ra + ^Iraim+l)' (87) 
{,& m=i,l,... 

Using similar techniques for the other cases, we find for the effective potentials 

V^(B = 0) = _l8 + 3C«(-3)-3Cfl(-l), (88) 

9 /77 5 /T7T , 11 11 ,11 

2V6- 2^ 10 +47-T^ + T 



V (1) (Sph) = -i-3v / 2-12v / 3 + |v / 6-|v / 10 + 77- ^7 + ¥ lo g( 2 ) 



+T l0 S (277) + Cf(-|, -3) + Cf(-|, -2) + GK-f, 1) 

+3 Cf(-|, -3) + Cf(~1, -2) - 5 Cf(-|, 1) (89) 

The pole term ^ has to be absorbed through a renormalization of the coupling constant. 
Adding the classical potential at the sphaleron we get 

V„fc(Sph) = — k 1 h finite renormalization. (90) 

env gl 2 4e v ' 

So the infinite part of the renormalization is 

1 1 11 

7 R ~7o + 12^ £ 



2 " .2 + (9!) 
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This is the standard result for SU(2) gauge theory in 1 + 3 dimensions as of course it 
should be: the renormalization is an ultra-violet effect and does not depend on the global 
properties of the space on which the theory is defined. 



3.3.4 General background field 

For a configuration along the tunneling path, c" = — uSf, d = 0, we make an expansion of 
the eigenvalues of the operators in u. The operators are: 

-AH AH = 4L 2 + 4uL ■ f + 2m 2 , (92) 
-A(w)(l -iV)A(w) = 4L 2 + 4uL ■ f + 2u 2 

-\u 2 (2 - (L ■ ff - L ■ f) Pf, (93) 

W(u) = 2L 2 + 2K 2 + 4uL-f+(Qu-2u 2 )S-f + 2u 2 . (94) 

From eq. (|9~2^) it directly follows that the eigenfunctions of the scalar operator A = 
—D i (B)D i (B) are the \lrriR; j rrij) with the eigenvalues 41(1 + 1) + 2u(j(j + 1) — /(/ + 
1) — 2) + 2m 2 . This allows us to write in eq. (|83|) 

oo l+l 

tr(e- sA ) = ]T (2l + l)e- s41 ^ ]T (2j + l) e -'^^-+D-'(W)-2)+^) ) (Q5) 

and to expand the second exponential in u. After performing the summation over j and 
the integration over s, the /-summations can again be expressed in terms of (^-functions. 



Most of the algebra was done using FORM [24]. The other scalar operator can be treated 



similarly. The eigenvalues can be written down from eq. (^3|) as well: they differ from the 
ones above only for 1 = 1. 

For the vector operator the situation is more complicated since the eigenvalues are not 
readily available. We write 

W(u) = 2L 2 + 2K 2 + W(u), (96) 



where the precise form of W(u) can be read off from eq. (|94j ). Since W(0) = we can 
treat W(u) as a perturbation on the operator 2L 2 + 2K 2 . As W commutes with the 
{L 2 , L 2z , J 2 , J z }, the basis to use is \ (ll)k,j). The dimension for the subspace in which we 
have to find the eigenvalues of W is given by the number of possible k values. Since this 
dimension is maximally three, it is of course possible to obtain the eigenvalues of W(u) 
exactly, but as explained earlier, it is more convenient to expand in u. The unperturbed 
eigenvalues within each sector are 21(1 + 1) + 2k(k + 1), so we need non-degenerate per- 
turbation theory to obtain the expansion in u for the eigenvalues X(j, k,l). A method that 



implements this particularly nice is due to Bloch p5 



Using MATHEMATICA |26| we computed the expansions for the eigenvalues. For this 
we needed the matrix elements of the operators L ■ T and S ■ T. To calculate these matrix 
elements, we need the 6-j symbols [27]. For three angular momenta t, s and I we have (cf. 
section |]) 

\t,(sl)k,j) = Y(-l) t+S+l+ n(^ + mk + l)]H t l S . q \\(ts)q,l,j). (97) 
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This gives 

(t, (Sl)k', j\S-f\t, (Sl)k,j) = -§( S (s + 1) + t{t + l))S k r k 

+lY / q( q + l)(2 q + l)[(2k + l)(2k' + ltA * * M j * * J- (98) 

The matrix elements for L-f follow from this and the identity J 2 = K 2 + T 2 + 2L-T + 2S-T. 
As in the case of the scalar operators, we expand the exponential in u and perform the 
s-integration. The summations over / can be expressed in ^-functions. During the various 
stages of these calculations care has to be taken for low values ofl,k and j: for certain values 
the Bloch result for the general eigenvalue X(j, k, I) will not be valid, since the subspace 
in which to perform the diagonalization may have less dimensions than the bulk value 
3 — \l — j\: some of the k values would be negative. Moreover, in the I = 1, j — 1 sector, 
which is a priori three dimensional, the k = modes have to be thrown out explicitly, 
since they are the d modes. Exclusion of the c modes is simply achieved by starting the 
summation over / with I = \- The result for the effective potential up to fourth order in u 
is 

V (1) (n) = -18 + 3 Ch(-3) - 3 Ch(-1) 

+ W (-H + h-T2^-Tk^+ l i log(2) + f log (^) 

-f c(§, -i) - \ C(| -i) + f C(I, -i) + 14 C(f , -i) 

+TU3)-TC«(5)-!C*(7)) 
+ u " (l§ + I ^ + 13 C(|, -1) - 9C(f , -1) - 106C(|, -1) 
-84 C(|, -1) + 2 C R (-1) - 23 £i(3) + 19 Cr(5) + 2 C«(7)) 
, 4 /863 , 1475 /K 53^/3 -,\ , 17>/5 i\ , 261>/7 i\ 



+47C(|, -1) - 66 C(f , -1) + f Ck(3) - 79 C«(5) 
+ ¥Ci?(7) + 5C«(9)). (99) 



V768 ^ 2048 v 4 S ^2' V ^ 4 

;|,-i)-66C(f,-i) + f 

z, 

4 

Here we eliminated the term in favour of the classical potential V^j defined by 

V C 1 = ^V d , (100) 

so 

V d ( M ) = \u 2 (2-u) 2 . (101) 

One checks that the pole term is absorbed by the coupling constant renormalization of 
eq. ©. 

Results up to tenth order in u were calculated. As can be seen from fig. |||, the expansions 
do not converge to the exact result at u — 1. This should come as no big surprise, since 
we have no reason to expect the radius of convergence of the expansion to be as large as 
one. Judging from the picture, one would estimate a radius of convergence of roughly 0.6. 
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To find the effective potential for larger u, we make a similar expansion of the effective 
potential around the sphaleron: u = 1 + a. The expansions of the determinants of the 
scalar operators can be obtained exactly as described above. The final /-summations can 
be expressed in terms of b) with b = —3, —2, 1. For the vector operator, we could in 
principle again use the Bloch perturbation method to obtain the eigenvalues A(j, k,l)(a), 
although the levels are now degenerate at a = . There is however an easier way. Suppose 
we have A = F + A with F = A(0) such that [F, A] = 0. This allows us to substitute in 
eq. (H) 



tr 



(e~ sA ) = tr (e- sF e~ sA ) = £ e~ sFi £ ^f- tr< (A n ) . (102) 



n=0 U - 



The sum over % is a sum over the eigenspaces of F, Fi is the corresponding eigenvalue and 
trj () denotes a trace within the eigenspace. For the operator W 

W(a) = 2L 2 + 2J 2 -2 + a(4 + 4L-f + 2S-T) + a 2 (2 - 2S • f ), (103) 

we indeed have that F = W(0) = 2L 2 + 2 J 2 — 2 commutes with W and we obtain 

lndet(W) = -^(^) £ £(2/+l)(2j + l)x 

/ ds s'l-i e -sm+i)+m+D-D J2 try . (104) 

n=0 n - 

Note that the sum over different fc values is now absorbed in the trace in the space. 
When taking this trace one has to deal with the same subtleties connected with low values 
of I and j as described before. The resulting effective potential is an expansion in a. Its 
behaviour is similar to that of the previous expansion: we have a radius of convergence of 
roughly 0.4 around the sphaleron u — 1. Using the fourth order expansion in u and the 
first order expansion in a (i.e. the value and the slope of the potential at the sphaleron), we 
can construct a polynomial in u of degree six that is a good approximation to the effective 
potential. 

Regarding the issue of the radius of convergence, one might think that also integrating 
out the (c, d) modes that are orthogonal to the u mode might result in better convergence. 
Also one expects the u — > 2 — u symmetry to be restored in this case. With the techniques 
described above, performing this calculation is straightforward, provided one remembers to 
properly adjust the gauge fixing procedure. The new expansions in u and a have roughly 
the same convergence behaviour as before, so there is no improvement on this point. The 
symmetry is however manifestly restored: the expansion in a contains only even powers of 
a. 

For the case of the vacuum and the sphaleron configuration, we could calculate the 
spectra of the operators exactly. For a configuration along the tunnelling path, the op- 
erators have less symmetries, and the evaluation of the spectra becomes harder. For 
B = B(c,d = 0) and especially B = B(c,d) we have even fewer operators that commute 
with the operators whose determinants we want to calculate. Nevertheless, we developped 
methods |18, ^] to calculate the functional determinants up to fourth order in c and d. 



These methods are based on eq. ( |102j ) and a generalization thereof. We postpone writing 
down the result of this calculation, which is the effective potential, until we have performed 
the renormalization. 
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3.4 The effective kinetic term 



To obtain the one- loop contribution to the operator B 2 , we perform the usual expansion of 
the path integral in Feynman diagrams. The subtleties related to the summation over the 
space-momenta were dealt with in the previous section. The diagrams needed are simple 
loop diagrams with one and two insertions of the operator A defined above. The particle 
in the loop is respectively a ghost, a scalar Q or a vector particle. The B 2 term comes 
from the diagrams with two insertions. 

As an example we will treat the diagram d 2 with two ghost propagators and two inser- 
tions of the operator A. The fermionic path integral is given by 



D'lp D'lp exp 



= J DtyD'ipexp — J tr 



9o 
1 

9o 



\fti (${-D,(B)(l-P)D„(B)} 
-d 2 - d 2 + 4L 2 } i) + $A 



The propagator in momentum space is given by 
-1 



(105) 



(106) 



k 2 + k 2 + 41(1 + 1) ' 

We will suppress the manipulations with the integrations over k e and x £ . The dimensional 
regularization can be obtained by adding k 2 to the denominators of all propagators and 
performing the integration 

d% 

WY 

Suppressing mi, rriR and m t , we find for the diagram with two insertions 



f d £ k £ 
J (27rK 



(107) 



d 2 = 



dkn dk n 
-1 



i v 



EE / dtdt ' 



k 2 +4l(l + l) 



(l\A(t')\l') e 



-iko(t-t') 



X 



k' l n + 41' (I' + 1) 



(l'\A(t)\l) e - ife o(*'-*). 



(108) 



x 



Noting that [L 2 , A(t)] = and substituting k Q = p, k' Q = p + q we write 
dpdq ^[dtdt'e-^ 
1 

p 2 + 41(1 + 1) {p + q) 2 + 41(1 + 1)' 
Using contour integration, one proves 

dp 1 1 1 



^ -ij^^jdtdt'e-^Hr^A^t' 

1 



/ 



2tt p 2 + C (p + q) 2 + C VC(4C + q 2 )' 
After performing the p integration and expanding in q 2 we obtain 
dq 

' i 



(109) 



(110) 



d 2 = -U — J2 [ dtdt'e-^'hiJA^A^' 



x 



E 



(-IT (q 2 ) 



2\n 



4(41(1 + 1))! t^o 4 " (4^ + 1))" 
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111) 



Next we substitute 
d 2 



Q 



dtdt' 



(112) 



and move these time derivatives to A by partial integration. The q and the t' integration 
thus become trivial, and we obtain 



d 2 = -| dtJ2Yl 

J l n=0 

„ oo 



.iy» 



1 



4« 
-1) 



(4Z(Z + l)) n+ ^ 

0n+g( e ) 



tr, (i (n) (t)i (n) (^)) 



4 7, 



tr, 



(iW(t)iW(t) 



(113) 

i n=o - W - 1))" - ' v 

The arrow indicates that we have performed the e-dimensional integrations that were 
suppressed in the notation. We introduced 



9n{e) 



d e x f 



^yr(n-|) 



(2tt) £ (1 + kf) n W*J T(n) ' ^ 14 ^ 

The n = contribution in eq. (|113|) can be related to the s 2 term in eq. ( 102 ), whereas the 
1 contribution gives us the B 2 term that we were looking for: 



n 



-i) + C(- 3 



The bosonic path integral is given by 

Z : 



(115) 



% 

D'Q D'Qi exp 



/ D'Qfj, exp \ [ tr(-2(D^)(fl)Q„ + Q M W^(fl)Q,,) 

\ [tr (Q {-d 2 -d 2 + AL 2 }Q 
+Qi {-0 2 -Bl + 2L 2 + 2K 2 }^ Q 3 + QoAQq + Q i W ij Q j 

-4Q ad{B j )Q j -2J Q -2JiQ i ) . (116) 

with J v = {D^F^ V ){B) (see appendix |BJ). The diagrams where the particle in the loop is a 
Qo particle have the same structure as those with a ghost particle. Apart from extra terms 
coming from the projector 1 — P in the ghost case, the contribution is precisely — ~ times 
the fermionic contribution. The diagram with a vector particle gives rise to some extra 
technical problems, but it can also be dealt with. There is also a diagram with an explicit 
dependence on B: it is the diagram with two insertions of the operator W j oc ad (-&;), one 
Qo and one Qi propagator. 

Adding up the different contributions, we obtain the one-loop contribution to the kinetic 
term c"c" + dfdf in the lagrangian: 



£«(c,d) = (c»c» + «) (t 



1247 



+ 



V2 + 



11 



11 



log(2) 



152 1 16 v " 1 12 e 24 7 + 12 

+T2 M^) + \ C(f , -i) + ii C(f , -i) + f C(|, -i) 

— | Cr(3) — J2 Cr(5) — | Cf(|> — ly 



117) 



The divergency will be absorbed by the same coupling constant renormalization of eq. (191 
that also made the potential part finite. 
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K — 


0.0566264741439181 


Ki = 


-0.2453459985179565 


K 2 = 


3.66869179814223 


K 3 = 


0.500703203096610 


K4 = 


-0.839359633413003 


K 5 = 


-0.849965412245339 


Kq = 


-0.06550330854836428 


K7 = 


-0.3617122159967145 


K$ = 


-2.295356861354712 



Table 1: The numerical values of the coefficients. 



3.5 Renormalized results 

We will use a renormalization scheme such that the renormalized kinetic part looks just 
like the classical term: 



eff = S + «) • 



29% 



1181 



Using eq. ( 117 ) this gives us the finite part of the renormalization (pT|): 
1 1 11 



gl + 12vr% 



11 1 ( L 



12tt 



(119) 



where k can be found in table |I[ This renormalization scheme can easily be related to 
other schemes like the MS or MS scheme. 

The finite, renormalised effective potential becomes 



2tt 



V e ff=^clM) + ^ 



with 



9r 



(120) 



V^(c) + V^(d) + K 7 tr(X) tr(y) + K 8 tr{XY), 
Ki tr(X) + k 2 det(c) + k 3 tr 2 (X) + k a tr(X 2 ) 
+k 5 det(c) tr(X) + /t 6 tr 3 (X), 



1211 



with the numerical values for Ki in table |l[ Note that the u 5 term in the effective potential 
along the tunnelling path uniquely determines the coefficient of the tr(X) det(c) term. The 
m 6 term can be obtained from combinations of the three independent invariants tr 3 (X), 
tr(X) tr(X 2 ) and tr(X 3 ). We choose to replace the w 6 term by tr 3 (X), which is the 
simplest of these from the viewpoint of the variational calculation. Note that not all of 
these coefficients are small, which means that the one-loop correction to the spectrum, to 
be calculated in section f|, may be substantial. 
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4 Calculating the spectrum 



4.1 Introduction 

In this section we will approximate the spectrum of the effective hamiltonian by applying 
the Rayleigh-Ritz method [213]. This is a variational method that consists of truncating 



some suitably chosen basis of functions and calculating the matrix of the hamiltonian 7i 
with respect to this basis. The results of the numerical diagonalization are the energy 
levels in various sectors. Using these, we find the excitation energies which correspond to 
the masses of the low-lying glueball states. We also estimate the window of validity of our 
effective model. 

Consider a hamiltonian 7i whose spectrum u{7i) is {jiti, {12, ■ ■ •}■ The Rayleigh-Ritz 
method gives upper bounds jl n > fi n . We will use the generalized Temple inequality |29|, |3(J 



to arrive at lower bounds for the levels. Let ip be a normalized trial wave function and 
define 

7 = (VW>, v = mn- 1 ) 2 \4>) = (ij\H 2 \*P)-j 2 - (122) 

Note that rj > and that rj = implies that ip is a true eigenfunction of Ti. With the help 
of the values of 77 we can derive lower bounds. We assume our variational basis to be so 
accurate that ji n +i > (p* n +i + An)/2. Under this and some other mild assumptions we can 



conclude 29, BQ 



2r7 . N 

[In - ~ — < /in < /V (123) 

We can treat the lowest order and the one-loop case on the same footing by writing 
f ( d 2 d 2 \ q 2 

n = -2{o^ + a^) +VM - s= h> (124) 

V(c,d) = V(c) + V(d) + \ 7 r 2 c r 2 d + X 8 mv 2 (c,d), (125) 
V(c) = A 1 r 6 2 + A 2 inv 3 (c) + A3r 6 4 + A4inv4(c) + A 5 r 6 2 inv 3 (c) + A 6 r c 6 . (126) 

The coefficients A, (which are functions of / and Kj) as well as the definitions of the used 
invariants can be found in table 0. We have also given the definitions of the invariants that 
occur in V 2 . For the lowest order effective hamiltonian, we just set Ki = 0, whereas for the 
one-loop effective hamiltonian we need to use the values for of table [1|. A sketch of our 
methods and some results for the truncated model were published before ||19|| . 



4.2 The variational basis 

We decided to use functions of the form 4>(r c , rd)Y(c, d) and to incorporate the boundary 
conditions in the (r c , r^) plane. Apart from the boundary conditions, we must also re- 
spect the symmetries of the hamiltonian as much as possible to obtain an optimal block 
diagonalization. Most of these symmetries, including the residual gauge symmetry, will 
be incorporated in the functions Y(c, d). After this, we will use the exact solution of the 
strong coupling limit of our hamiltonian problem as a guide to obtain a useful set of radial 
functions. 
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Ai 




A 2 


-7 + 3^2 


A 3 


= «3 + K A 


A 4 


= j — 8K4 


A 5 


= |«5 


A 6 


= K 6 


A 7 


= K 7 + Kg 


A 8 


= f - 6k 8 



inv 3 
inv 4 
inv 2 


[c) 
[c) 

'c, d) 


= 3det(c) 

= |(tr 2 (X)-tr(X 2 )) 

= |(tr(X)tr(F)-tr(XF)) 


inv 6 


[c) 


= mv s( c ) 


inv 7 


[c) 


= inv 3 (c) inv 4 (c) 


inv 8 


[c] 


= mv 2 (c) 



Table 2: The coefficients for the hamiltonian and the definition of several invariants. 



4.2.1 Symmetries 

The hamiltonian T~C(c, d) is invariant under the transformation c — > ScRi, d — > SdR 2 with 
S,Ri,R2 G SO (3) and under the interchange c <->• d. The generators of left- and right 
multiplication are L^, Lf , and L|. These are su(2) angular momentum operators and 
we have for instance 

d 



( L c)=-^-^. (127) 

The following set of operators commutes: 

{W,J s ,J^,(Lf) 2 + (^) 2 ,P}. (128) 

V is defined by Vf(c,d) = f(d,c). On S 3 it corresponds to the parity (n ,n) <-> (no, — n) 
(cf. eq. (^)). The operator J 5 = Lf + L| implements constant gauge transformations: we 
have to demand (J 5 ) 2 = for physical wave functions. The operator J R = Lf + is the 
rotation operator. The spatial symmetry group is SO (4) and these symmetries cannot be 
simply divided in independent sets of translations and rotations. The operator on S 3 that 
corresponds to J R is L\ + L 2 , and the corresponding symmetry on SU(2) (cf. appendix [Xp 
is g ^ g\ X gg\- Writing g = n ■ a and g\ — n' ■ a we have 

n-a — ► n^(n' ■ a)a^(n' ■ a) = n cro — riiV- 1 (n')aj. (129) 

The symmetry on S 3 is thus seen to leave Uq invariant and to perform an SO (3) rotation 
on n. It is hence a rotation for points around the north and south pole of S 3 . The 
operators L\ and L 2 do not leave any point invariant and cannot be interpreted as rotation 
operators. Only their sum has this property and the different sectors (j = 0, 1, . . .) under 
the symmetry J R correspond to scalar glueballs, vector glueballs, etc. 

To prepare the ground for the block diagonalization of the hamiltonian, we divide the 
function space in sectors characterized by the quantum numbers j, m, h(h + 1) + h(h + 1) 
and p corresponding to the operators of eq. (128). Note that for low values of l\ and 



l 2 there is a one-to-one correspondence between unordered pairs (h,h) and the numbers 
l\{lx + 1) + l 2 {h + 1)- Since the spectrum of Ti is independent of the azimuthal quantum 
number m, we use the notation lil 2 j-eveia and /i/ 2 J-odd to denote the various sectors. The 
lowest energy level in the scalar (j = 0) sector is the vacuum. Energy differences with 
respect to this level will be the glueball masses. 
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Operator 


(L S C ) 2 










Eigenvalue 


1.(1, + 1) 


m s 


l r (l r + 1) 


m r 


—L(L + 7) 



Table 3: Behaviour of the functions (c\L; l s , l r , r; m s , m r ). 



4.2.2 The angular sector 

We begin with constructing functions of c that are eigenfunctions of the following set of 
commuting operators 

{ A* (130) 

Note that the space of the c makes up an S 8 . From appendix [A] we know that the eigen- 
functions of the spherical laplacian will be the homogeneous harmonic polynomials in 
c. An orthonormal basis of functions of c is given by the set {(c\L; l s , l r , r; m s , m r )}. Each 
of these functions is a homogeneous harmonic polynomial of degree L in c. Its eigenvalues 
under the various symmetries are collected in table [| We used the operators L R and 
to further classify these functions, but for higher values of L we need the extra label r for 
the remaining degeneracy. This degeneracy is absent for low values of L, shows up first for 
(L, l s , l r ) = (4, 2, 2) and proliferates for high values of L. 

Explicitly, we have that (c|0; 0, 0, 1; 0, 0) is just the constant function. For L = 1 we 
have the nine functions c". With the help of the raising and lowering operators L± and 
L± we arrive at the nine functions (c|l; 1, 1, l;m s ,m r ) which are the appropriate complex 
linear combinations {c+, c° , c^, Cq, ■ ■ ■ ,cz} of the variables c\. We have e.g. £fc+ = — c+ 
and Lfc+ = cfj\ 

To construct all representations (L,l s ,l r ,r) for given (L,l s ,l r ), we proceed as follows. 
We construct the eigenspace Af of the operators L§ and Lf. This is the space span by 
all monomials of degree L in the complex c variables that have the eigenvalues l s and 
l r under the operators L§ and Lf. Within Af, but now regarded as functions of c in 
stead of c, we construct the intersection of the kernels of the operators L_L+, L^L+ and 
(cfc"A). The first two operators leave us with only those combinations of the monomials 
that have the prescribed (L s ) 2 and (L R ) 2 eigenvalues; the laplacian imposes the constraint 
that the polynomials must be harmonic: A = 0. If this intersection consists of more than 
one function, we use explicit Gram-Schmidt orthogonalization to arrive at the functions 
(c\L;l s ,l r ,r;l s ,l r ). With the lowering operators, we trivially construct the rest of the 
functions (c\L; l s , l r , r; m s , m r ). 

Most calculations were done in MATHEMATICA, but we used C to perform the con- 
struction of the intersection of the kernels of L^L+ and L^L^. For this we used the explicit 
projectors on the kernel 

p -°%i h +i)Vi) L -4 (i3i) 

and the analogous expression for P t . Since all calculations had to be done exactly, we used 
arithmetic with (large) integers. In this way we explicitly constructed the representations 
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Operator 


(^f) 2 


(4 5 ) 2 


(J s ) 2 




Ac 


Eigenvalue 


Isih + 1) 


1.(1. + 1) 








-Li(Li + 7) 




Operator 






(J R ) 2 


J£ 


A rf - 


Eigenvalue 








m 


-L 2 (L 2 + 7) 



Table 4: Behaviour of the functions (cd\j, m, l s ; L±, h, t\\ L 2 , k, t%). 

(L,l s ,l r ,r) for L < 10 and l r < 2. As will be apparent from the sequel, we did not need 
higher values for l r . 

We construct functions of both c and d by using the familiar rules of adding angular mo- 
menta. Let i denote a representation (L;l s ,l r ,T) and consider the functions (c\ii,m s ,m r ) 
and (d\i2', Tn' s , m' r ). Using Clebsch-Gordan coefficients, we can define a function Y ni2 (c, d) 
which is an eigenfunction of J R and of J s . We will limit the construction to functions with 
J s = 0, as required by residual gauge symmetry. This implies that the functions of c and 
d need to have the same l s , which restricts the possible combinations of i\ and % 2 . The 
resulting function is given by 



Y llt2 (c, d) = (cd\j,m,l s ;L 1 ,l 1 ,T 1 ;L 2 ,l2,T 2 ) 



Is h I 



■2 



h k i 



x 



= E E E (-i) ,1 - h+B v j y + 1 

ms= _ lsmi= _ hm2= -i 2 \m 1 m 2 -m 

, (c\Li, l s , li, n; m s , mi) (d\L 2 ; l s , k, t 2 \ -m s , m 2 ). (132) 

y ll a + 1 

Its eigenvalues under the various symmetries are collected in table Note the behaviour 
under parity: 

-p^ii«2 _ t ^ii+i2+iy«2w (133) 



4.2.3 The radial sector 

The strong coupling limit of our hamiltonian problem consists of the eigenvalue problem for 
the kinetic part of the hamiltonian of eq. (fl9|). If we assume a solution (^(r c )(^(r^)y* li2 (c, d), 
the reduced one-dimensional eigenvalue problem becomes 

(&_ + x _ (12 + L(L + 7)) 
\ dr 2 r 



2 + A - ^ A >± v{r ) = 0, (134) 



whose regular solution is tp^\r) = 7^J3+L(7 r ), with j p (z) the spherical Bessel function of 
order p and A = 7 2 . The eigenfunctions of the kinetic part of the hamiltonian are thus 
given by 

^\r c )^\r d )Y^(cJ), (135) 
with the energies given by E = + 7I). 
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During the variational stage of the calculation, however, the use of spherical Bessel 
functions of different order will lead to a large number of integrals. Therefore we take the 
radial functions to be independent of Li and L 2 and define 

^%(c, d) = ^ 71 (r c )^ 72 (r d )y ill2 (c, d), (136) 

with y? 7 (r) = ^rj^r). These functions are not eigenfunctions of the kinetic part of 
the hamiltonian and they will have discontinuities of the following kind. For r j 0, we 
have that y? 7 (r) ~ r 4 . The wave function \l/ behaves as r~ 4 ip 1 (r)(c\L; l s , l r , r; m s , m r ), 
and this function will be discontinuous at r = 0. The variational functions thus have 
discontinuities at r c = and = 0. These discontinuities form a set of measure zero, and 
a variational calculation will not feel them. Since the functions fy(r) are the eigenfunctions 
of the reduced one-dimensional problem with L = 0, they still constitute a complete set of 
functions. 

The behaviour under parity can be obtained from eq. ( |133j ) and is given by 



v^ 2 = ^) h+h+3 ^- ( 137 ) 

Taking even and odd combinations gives 

V#£(c,d) = %%{c,d) + p ^%(d,c). (138) 

We implement the boundary conditions eq. fl5?D and fl58|) for 9 = by imposing the 
following conditions on 71 and 72 : 

p=-l : Vll (^>) = f y A^)=«, (139) 

p=1 : ^^ ( V3) = *Z^) (V 3) = 0. (140) 

These conditions are expected to be accurate as long as the wave function transverse to the 
sphaleron path (near the sphalerons) is predominantly in its ground state. The case 9 = tt 
can be treated along the same lines as 9 = by interchanging the boundary conditions for 
the cases p = 1 and p = — 1 . 

The exact strong coupling results for the lowest levels in some sectors for 9 = are 
collected in table [5|. These values are obtained by imposing the boundary conditions above 
on the true eigenfunctions of eq. the values of 71 and 72 are hence dependent on Li 



and L 2 respectively. 

For general 9 we multiply (c, d) with a phase factor exp(i9a(r c , r d )). The func- 
tion a is a kind of Chern-Simons functional that gives the correct behaviour to the wave 
function under large gauge transformations. The resulting functions no longer have well- 
defined parity, but they do obey the general boundary conditions for suitable a. Also 
the hermiticity of the hamiltonian for these functions can be checked explicitly. Sufficient 
conditions on a are: a{r C) rd) = —a(rd,r c ) and a(v^3, 0) = |. We choose 




(141) 



For (3 — > 00 we approach the situation that the phase factor over the entire edge is constant 
and equal to e %e . But already for the choice (3 = 2, the boundary conditions at the 
sphalerons are taken into account properly. 
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sector 


\j,m,l s ;L 1 ,l 1 ,Ti;L 2 ,l2,T2) 


7i 


72 


E 


000-even 


|0, 0,0; 0,0,1; 0,0,1) 


1.9786 


1.9786 


3.9149/ 


000-even 


|0, 0,0; 0,0,1; 3, 0,1) 


1.9786 


4.1215 


10.4508/ 


000-odd 


|0, 0,0; 0,0,1; 0,0,1) 


4.0345 


6.0143 


26.2246/ 


112-even 


|2,m,l; 1,1,1; 1,1,1) 


2.7406 


2.7406 


7.5108/ 


112-odd 


|2,m,l;l,l,l;2,l,l) 


4.7242 


5.4016 


25.7476/ 


022-even 


|2,m,0;0,0,l;2,2,l) 


1.9786 


3.4456 


7.8936/ 


022-odd 


|2,m,0;0,0, 1;2,2, 1) 


4.0345 


5.4016 


22.7272/ 



Table 5: Strong coupling limit: lowest energy levels in some sectors. 
4.3 Matrix elements 

In this section we will describe the calculation of the matrix elements of TL and of Ti 2 . 
The hamiltonian Ti = K, + V is given by eq. (fl9|) . We let \n) denote the normalized basis 
function: 

\ n ) K te 2 2 M) +p(-l) h+l2+j ^(c,d))exp(iea(r c ,r d )), (142) 



with ip l ^ 2 (c,d) given by eq. (|136|) . We take a from eq. ( |141| ) with ft = 2: a(r c , r d ) = 

Consider first the matrix element (n'\H\n). For the potential energy V the phase factor 
exp(i9a(r c , r d )) cancels against its complex conjugate. For the kinetic operator we obtain 

exp(-i6a(r c , r d )) K (V^fo d) exp(i6a(r c , r d )) 

f ( 2 ^ 2 £i(£i + 7) , L 2 (L 2 + 7) 

77 7i + 7 2 + 

2 \ r 



c 



+ 2 ^ + ^WO^fa)) ^"(c, 4 (143) 

To apply Temple's inequality, we need the matrix elements of Ti. 2 . Using the hermiticity 
of Ti. we write 

(n'\H 2 \n) = ({JC + V)n'\{JC + V)n) 

= -^{K,n'\K,n) + (n'\VKL\n) + ^(n'\V 2 \n) + hermitian conjugate. (144) 

From these expressions we can read off which matrix elements we have to calculate. 

Since the potential V(c, d) is invariant under and as well as under J 5 , all the 
terms in V and V 2 are of the form 

A = (cd\0, 0, l s ; Li, 0, n; L 2 , 0, f 2 ). (145) 
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By virtue of the construction in eq. (|132|) all the angular integrations over c and d can 
be reduced to angular integrations in H 9 . These integrals can be expressed in terms of 
reduced matrix elements. The details on this procedure can be found in appendix [C[ 
The radial integrals we need to compute are 

.7(71,7,7) = / cirr n <£y(r)v3 7 (r), (146) 



JiM,i) = jf 3 drr n <p Y (r) f r ^^ 7 (r)j , (147) 



Ju{ n il'il) 



J Q Z drr n ^0y(r)j [r^{r)\ , (148) 



where (p^ is normalized to one, and 

Vl ( r ) = jrj 3 (jr) = f(jr), (149) 
f(z) = (l-^)cos(z)+(-^ + p)sin(^). (150) 

Consider the integral for non-normalized functions and let a = y/S'j and a' = v^3~7'- For 
an explicit value of n, but unspecified values of 7' and 7, we write 

J(n) = [ 3 drr n ip Y (r)^(r) = (V3) n+1 C dx x 11 f (a' x) f (ax) . (151) 



We express products of sines and cosines of ax and a'x in sines and cosines of sx and vx 

with s = a' + a and v = a' — a. The integrand as it stands is regular at x — 0, but 

individual terms need not be. We therefore introduce a cut-off £, after which the integral 

is a sum of functions 

1 ,1 
dxx m cos(ax), sin(m, a) = / dxx m sin(ax). (152) 

Using recursion formulae obtained from partial integration allows us to eliminate the func- 
tions cos(m, a) and sin(m, a). After taking the limit e j 0, the result is an expression in 
the functions sin, cos, Si, Ci and In of the variables s and v. By substituting values for s 
and v and performing the normalization we construct tables for J(n, 7', 7). Note that we 
have to be careful for a' = a: we have to take the limit v — > before substituting actual 
values. 

The other integrals Ji(n) and Ju(n) can be obtained as follows: 

rs/3 ( d \ 

Jl{n) = drr n ip Y (r) I r— ip y (r) I 

= (y/3) n+1 J 1 dxx n f(a'x) (x^-f(ax)\ 



J(n). (153) 



s — v ( d d 
2 \ds dv 
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4.4 Weak coupling expansion 

In this section we will perform the perturbative calculation of the energy levels of the 
hamiltonian. The results will serve as a check on the variational calculation in the regime 
of small coupling constant. Starting from the strong coupling basis, we could also develop 
a perturbation theory in 4, but there are two reasons why this is less interesting. First, 
we are interested in the region where we just start to see the deviations from perturbative 
behaviour in /. Second, since our variational basis is in essence a strong coupling basis, the 
reproduction of the strong coupling limit does not give a strong check on the variational 
calculation. 

For the weak coupling limit, we rescale x c = \Jjr c an d x& = Jj^d and we split the 
hamiltonian of eq. (|49| ) in H = Ho + Hi, with 

( d 2 d 2 (I 1 \ 1 1 \ 2 2 

H ° = -(^ + ^- 12 ^ + ^j + ^ + ^l + ^ + ^ (154) 

fti(c,d) = Hi(c)+Hi(d) + X 7 U\ x 2 c x 2 d + X 8 U) x 2 c x 2 d mv 2 (c,d), (155) 



Wi(c) = (|A 1 -ljx c 2 + X 2 x 3 c mY 3 (c) + (|) 2 A 3 x^ 




^ X 4 xfmv 4 (c)+ (^J A 5 ^inv 3 (c) + ^J X 6 x b c . (156) 
Diagonalizing Ho leads to the following one- dimensional eigenvalue problem: 

(* + A- 12 + L ( L + 7 > -^V M = 0. (157) 



dx 2 x 



y(x)=e~* x x L +\F x {-- -,L + %,x 2 ), (158) 



The regular solution of this equation is [|3T| 

.1-2 T.J.A ^ , X — 2L — 9 

1 

with iFi the confluent hypergeometric function. As in the strong coupling case, we can use 
the boundary conditions at x — to discretize the possible values for A and hence the 
energy. This would give us an alternative to the strong coupling basis constructed earlier. 
Disadvantages of this basis are that the functions themselves depend on the coupling 
constant /. Also the matrix elements of operators between two confluent hypergeometric 
functions are hard to calculate. For the variational method, we will therefore use the strong 
coupling basis. 

We can however use the weak coupling basis for perturbation purposes. For / j the 
location where the boundary conditions have to be imposed moves away to infinity, which is 
equivalent to saying that the wave function, in the original coordinates, is strongly localized 
around c = d = 0. We can replace the boundary conditions at x = ^Jj by demanding 
normalizability of the functions <p(x). This means that the iiq function must reduce to a 

polynomial, which is the case if its first argument is — n with n = 0, 1, The \Fi function 

then reduces to a Laguerre polynomial and we obtain the discretization A = 9 + 2L + An. 
The weak coupling eigenfunctions are thus 

« 2 = ^^^xD^^V^x^^x^^e-^e-^^-^c,^, (159) 
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and the energies are given by 

E = 18 + 2(L 1 + L 2 )+4(n 1 +n 2 ). (160) 

We use Blochs formulation |25], [l|] of perturbation theory to calculate the weak coupling 
expansion of the energy levels of the hamiltonian. For this we need the matrix elements of 
Hi w.r.t. the basis of eq. ( |159| ). The angular integrations are treated in appendix y. The 
radial integrals take the form 

r dxx n e- x \ L ' +A x L+i L^ +l \x 2 )L { n +l \x 2 ). (161) 
Jo 

We calculate these integrals by substituting the explicit form of the Laguerre polynomials. 

Using orthogonality properties of the Laguerre polynomials and of the angular func- 
tions, one can show that the number of intermediate levels that we have to take into 
account is finite. To be more precise: given the values E and L for an unperturbed level 
and the requested accuracy in terms of powers of the coupling constant, one can calculate 
maximum values E max and L max beyond which there will be no contribution up to the 
requested accuracy. 

We will give the perturbative results for those levels which will turn out to correspond 
to the lowest glueball masses: 

io , W9 . 9 \ , f 2 ( 63 9 2 3 

+ f k 3 + f «4 + § « 7 + § «s) + O (f) , (162) 

oo , * ( 9 , 11 \ , ,2 /903 , o 11 2 13 

22 + / (-g + T KiJ +/ Ui2 + 3 «i-l6«l W 2 

i 165 , 105 ,135 45 \ rn f f 3\ 

+ -8-/^3 + — «4 + l6" K7+ 16 Ks ) +°(f J' ( 163 ) 
00 , f ( 13 , 11 \ , f 2 /1063 ,13 11 2 13 

, 165 ,105 ,99 , 33 \ , /0 / ,3\ ( -,n A \ 
+ — K 3 + — K 4 + Jq k 7 + Kg) + £> (^/ J , (164) 

oo , r /13 , 11 \ , f 2 { 3167 ,25 11 2 7 

22 + f {24 + T^)+f (4608 + 24 ^ ~ 16 Kl " 8 K2 

+*f K 3 + f K4 + i| K7 + |I K8 ) + o(/3). (165) 



With the results of the previous sections, we are in a position to actually perform the 
variational calculation. A FORTRAN programme constructed the truncated matrix for the 
hamiltonian and performed the diagonalization using routines from the Nag-library and 
the Eispack library. The output consisted of the upper bounds jl n for the lowest energy 
eigenvalues and of the corresponding values for r\ n (used for computing lower bounds). 
We also produced graphical representations of the wave function. These can be used in 
verifying a posteriori whether the assumption is true that the boundary conditions are only 
felt at or near the sphalerons. 



^000-even 



Ei(l) 

^000-even 



^000-odd 



^112-even 



4.5 Results 
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4.5.1 Lowest-order hamiltonian 

We start by considering the lowest order results. After calculating the lowest levels in 
various sectors, we find that the lowest glueball masses are to be found in the sectors 
000-even, 000-odd and 112-even for respectively the scalar and the tensor glueballs. These 
levels, including the lower bounds and the perturbative expansions, are plotted in fig. |3|. 
Note that the crossing of levels only happens for levels that are in different sectors. States 
within the same sector would exhibit avoided level crossing due to level repulsion. 

For small coupling, the deviations of the wave function from the perturbative one will 
be small. To reproduce the energy eigenvalues correctly, we can suffice with a relatively 
low number of angular functions, but to be able to accommodate the strong localization 
around c = d = 0, we need a large number of radial functions. For larger values of / we 
can reduce the number of radial functions, but we must increase the number of angular 
functions. The data we used was obtained by combining different runs with various choices 
for the truncated basis as controlled by the parameters A^ rac j (number of radial functions) 
and L sum (limit on L\ + L 2 ). 

In the 000-even sector we maximally used N Ta( ^ = 150 and Lsum = 10 for small /, 
and iV" rac } = 50 and Lsum = 14 for larger values of /. These values correspond to bases 
consisting of over 3000 vectors. Increasing the number of basis vectors becomes quickly 
limited by the amount of free memory available in the computer. For the case 9^0, the 
hamiltonian is complex and we need a factor of two more memory. In the 112-even sector, 
the same number for Lsum resulted in higher numbers of angular functions. Here we went 
up to iV rac j = 150 and Lsum = 8 for small /, and N ra( ^ = 30 and Lsum = 12 for larger 
values of /. 

We did most of the calculations on Sun workstations. For some values of / we wanted 
more accurate results and used computer time on a Cray. We also developed a pruning 
technique: after diagonalization, we examine the variational wave function obtained for, 
say, the ground state. We replace those basis vectors (roughly ten to twenty percent) 
that have a small coefficient in this wave function by new basis vectors and repeat the 
diagonalization. Repeating this procedure allows us to sample a larger variational basis, 
but we do not have the guarantee that the results will improve in this way: it is possible 
that the true wave function has substantial overlap with very many vectors in our basis. 

Returning to our calculation, the lowest-lying scalar (j = 0) and tensor (J = 2) levels 
are found in respectively the sectors 000 and 112. For 9 = 0, the vacuum corresponds to 
the ground state of the 000-even sector. The scalar glueball + can be identified with the 
first excited state in the 000-even sector, the tensor glueball 2 + with the ground state in 
the 112-even sector. Note that at / = 0.6 these levels cross, thus making the scalar glueball 
heavier than the tensor. We do, however, not expect our model to be valid anymore for 
these values of the coupling constant. Taking differences of the energy levels with respect to 
the ground state in the 000-even sector gives us the masses of the low-lying glueball states 
(fig. |j). To check on the strong coupling results, we also performed variational calculations 
for large values of /. The comparison with the results of table [5] can be seen in fig. |5|. 

For small /, there is virtually no dependence of the masses on 9, whereas for larger 
/ the 9 dependence becomes bigger. This is shown in fig. |6|. For 9 = ir we can impose 
the boundary conditions exactly, that is, without using the trick with the Chern-Simons 
like operator, but with the same method that we imposed the 9 = conditions. The 
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variational results obtained in this way are also plotted in fig. [| Remember that a 9 
dependence is a sign that the boundary in the (r c ,r^) plane is felt. Our model is valid 
for values of the coupling constant at which only the boundary conditions at and near the 
sphalerons are felt. Checking this a posteriori with the help of plots of the wave function 
(fig. |7|) indicates that / should not be larger than roughly 0.5. These plots were obtained 
as follows. Consider the function 

|^(r c ,r d )| 2 = J dcdd \i/;(c,d)\ 2 , (166) 

which is a measure of the probability distribution in the (r c , r<j) plane. Dividing this 
\'4 ) ( r c, r d) \ by r^rj, we obtain a function with the expected behaviour of the true wave func- 
tion: it is localized at the sphalerons and decays exponentially in the transverse directions 
(see fig. 0). Note that the characteristic width grows with increasing /. Although the low- 
est barrier is the sphaleron, the measure factor r\r\ causes the configurations that are close 
to the sphaleron but have a somewhat higher energy to make the dominant contribution 
to the tunnelling. The relevant parameter here is the characteristic decay length of the 
wave function, which in turn is determined by the rise of the potential in the transverse 
directions. 

Although we are primarily interested in the case 9 = 0, an appreciable dependence on 
9 for a certain value of / shows that the non-perturbative influence of the boundary has 
become important. To explain the fact that the spectrum is not exactly periodic in 9, note 
that our implementation of the 9 dependence (eq. ( |141| )) only has this periodicity in the 
limit /3 — > oo. The volume effect described above implies that the relevant distribution 
iV'tjcj ?"d)| samples a piece of the boundary over which the phase difference already starts 
to depart from e td . This also implies that the results using the exact implementation of 
the boundary conditions for the case 9 = it can be expected to differ from those using the 
method of eq. ( |141| ). Increasing (3 would improve the periodicity properties, but would not 



result in a better effective model: although it would be less apparent, the wave function 
would still be spread out over the boundary. 

4.5.2 One-loop hamiltonian 

We present basically the same plots for the one-loop hamiltonian as for the lowest order 
hamiltonian. As before, the lowest glueball masses are to be found in the sectors 000-even, 
000-odd and 112-even. These levels, including the lower bounds and the perturbative 
expansions, are plotted in fig. |8]. Note the avoidance of level crossing between the ground 
state and the first excited stated in the 000-even sector, which gives rise (cf. eq. (|123|) ) to 
relative large errors for even the ground state. 

Especially in the 112-even sector, the values obtained for rj are rather large. We like to 
stress here that the lower bounds as obtained by Temple's inequality are rather conservative 
and that the actual error is often much smaller. This insight is gained in studying toy 
models, but the effect can also be observed in the case at hand. Going from L sum = 12 to 
Lsum = 14 in the 112-even sector, the upper bounds shifted very little, whereas a reduction 
of the 7] value by a factor of two was achieved. 

It is to be expected that for the one-loop hamiltonian a larger number of basis vectors is 
required, since the potential has a more complicated structure in this case. Also the values 
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of the coefficients in the potential may cause the large values of rj. When investigating 
the influence of Aj on the variational results, we find that the large value of A§ seems to 
cause the problems. This coefficient corresponds to the operator inv 2 (c,d) (see eq. (|125|) 
and table |2|), which is the only operator in the potential that can change the l s value of a 
basis function. Apparently, our basis could not be chosen so large as to yield high-precision 
results when the coefficient of this operator becomes large. 

Taking differences of the energy levels with respect to the ground state in the 000-even 
sector again gives us the masses of the low-lying glueball states: see fig. |9|. Fig. [K] shows 
the 9 dependence of the scalar glueball. In view of the remarks above, we left out the 
lower bounds in these cases. We can again study the localization properties around the 
sphaleron: see fig. |TT[ 



4.5.3 Discussion 



Using the Rayleigh-Ritz method we can determine the spectrum of the effective hamilto- 
nians. The use of Temple's inequality gives us confidence that our results are accurate, 
especially since experience tells us that the actual error is usually much smaller than the 
conservative estimates based on the values of rj. The results are also consistent with the 
strong coupling limit and with the weak coupling perturbative expansion. 

Let us focus first on small values of /. Here the boundary conditions are not felt yet, 
and the variational results are in accordance with the results from perturbation theory. 
When / grows, we see the effect of the boundary conditions in field space set in. For even 
larger /, the wave function has spread out over the entire (r c , r^) plane, and the model has 
lost its validity. It is however possible that this already happens before this point because 
of a break down of the adiabatic approximation. In order for the adiabatic approximation 
to be valid, the wave function for the transverse modes must still be in its ground state. 
To check this for the modes outside the (c, d) space is hard, but transverse to the tunneling 
path we can examine the transverse modes within the (c, d) space. Fig. ^ and |TT]do not only 
show the localization of the transverse wave function at the sphaleron, but also indicate 
that this transverse wave function is in its ground state. We therefore assume that the 
adiabatic approximation is valid, and that values of the coupling constants for which our 
effective model is useful are determined by the spreading out over the (r c , r^) plane. 

From these considerations, we derive the following windows in the coupling constant. 
For the lowest order hamiltonian, the interesting window is between / = 0.3 and / = 0.5. 
For the one-loop case, it is between / = 0.2 and / = 0.3. 

One of the issues raised above was the level of localization of the wave function around 
the sphaleron. This is related to the assumption that only the boundary conditions at and 
near the sphalerons are felt. We argued that this was determined by the rise of the potential 
in the transverse directions. The one-loop correction to the tr(Y) term in the potential at 
the c sphaleron, which can be expressed in Ki, k 7 and k 8 , is such that it results in a lesser 
degree of localization. In both the lowest order and the one-loop case, a strong localization 
of the wave function around the sphaleron is not realized and the boundary conditions 
are not felt exclusively at the sphalerons. This should not come as a big surprise since 
it is unnatural for an eigenvalue problem that the spectrum is determined by boundary 
conditions at a finite number of points. Typically the conditions at a boundary with 
codimension one determine the spectrum. This is illustrated by our problem: we feel the 
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conditions at a finite part of the boundary, determined by the localization potential in 
the transverse directions. As long as the coupling is low enough, our assumption that the 
boundary conditions are only felt close to the sphalerons is true and the results are valid. 



5 Conclusions 

The goal of this project was to calculate the glueball spectrum on S 3 . In particular, we 
set out to study the effects of the topology of the configuration space on the spectroscopy. 
These effects can be seen as the result of the presence of instantons. 

It is important to emphasize one should not expect our results for the spectrum to be 
accurate for large volumes. For large volumes the effects of the non-trivial topology and 
geometry (curvature of the configuration space, not to be confused with the curvature of 
S 3 ) become too strong to be described by the effective theory. Within the effective theory 
we clearly observe that at large coupling the wave functional will start to feel more of 
the boundary of the fundamental domain than just the neighbourhood of the sphaleron. 
However, it has been the main aim of this study to demonstrate that instanton effects on 
the low-lying spectrum are large, but calculable as long as energies remain close to the 
sphaleron energy, where nevertheless semiclassical techniques will completely fail. 

The fact that ttiq+R is decreasing, down from tuq+R = 4, is what should be expected for 
the following reason. A rough estimate for where one would expect instantons to become 
relevant, based on what one finds on T 3 , would be tuq+R ~ 1.4. Here we equate the largest 
geodesic distance on a torus of size L, \/3L/2, to its value on S 3 , ttR, or 

L = -j=R, (167) 
and we use that on T 3 instantons are relevant for z = m +L > 5 |32|, [33|. Furthermore, 



we assume that the scalar glueball masses are roughly equal in both geometries at these 
volumes. These low values of m +R are not reached in the tree level approach, but we 
do reach them in the one-loop case which shows that it was necessary to include the one- 
loop corrections. This regime of masses occurs for values of the coupling where we expect 
our model to still be valid. Specifically, m +R = 1.4 corresponds to / = 0.28. At larger 
couplings, we clearly see that the wave functional feels too much of the boundary of the 
fundamental domain. For / = 0.4 this is dramatically clear from fig. [XT], where the wave 
function is seen to probe unacceptable regions of the boundary of the fundamental domain 
to remain a good approximation to the full wave function. From fig. |] we see that the 
scalar and tensor mass even cross around / = 0.33, which is certainly unacceptable for the 
full theory. Clearly we have pushed the model passed its region of validity for / > 0.3. 

Of course, at some point m +R has to start to rise again, and when m + reaches its 
asymptotic infinite- volume value, tuq+R grows linearly with R. Both the truncated and the 
one- loop results show that m +R exhibits a minimum, after which it starts to rise again. 
It rises linearly in / for / — > oo, as follows from the strong coupling results in fig. [|, which 
are also valid for the one-loop case. This however does not mean that we are approaching 
the infinite- volume limit, because our effective model is certainly not fit to describe this 
regime. Moreover, it is clear that no statements can be made on the R dependence of / 
for these large couplings and volumes. 
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Other indications that our results are in the domain of expected validity are that at 
/ = 0.25 the tensor to scalar mass ratio is given by 1.5, rising from 1 at zero coupling to 1.8 
at / = 0.28, see fig. |9|. For a torus geometry one finds a similar range. Below z = m +L = 5 



the tensor is split into a doublet at ~ 0.9m + and a triplet at ~ 1.7m + |32], |34| (due to 
breaking of rotational to cubic invariance), whereas these states seem to merge at z > 7 
into a degenerate quintet with m 2 + = 1.5m +. 

We assume that in intermediate volumes the two-loop (3 function can be used to convert 
the dependence on the relevant coupling constants to the dependence on the radius R or 
the size L. Thus, for R we have (/ = g 2 /(27r 2 )) 



(A a fl) a = Wr" ^> A = ^?. ^ = ^> ( 168 ) 



- hi 11 & _ 102 

and the same formula holds for A^L in terms of the minimal subtraction coupling gus- If 
we, as usual, set the scale by a string tension of a = (425MeV) 2 and use that y/a/m + rs 0.3 
for T 3 , we have that in physical units z = m +L = 5 corresponds to L = 0.7 fm and hence, 
using eq. ( |167| ) , to R = 0.19 fm. As derived earlier, this corresponds to / = 0.28 on the 
three-sphere, and eq. ( |168f) gives A R = 1.3 fm -1 . On the torus, z = 5 corresponds to 
gMS — 2.42 0. Relating this to L = 0.7 fm gives Aj, = 0.39 fm -1 . In this way we can 
compare our result for m 2 + / m + as a function of the volume to the doublet E + and the 
triplet T 2 + on the torus, as is shown in fig. IT3. The vertical bar on the right indicates the 



range of lattice Monte Carlo values [[34]] for the E + and T 2 masses (equal within errors) 
at L = 1 to 1.5 fm. 

Another way to relate the length scales above to a value of the coupling constant uses 



the definition of the running coupling constant in [|35|] and the relation between this coupling 
constant and the one in the MS scheme. Proceeding in this way relates R = 0.19 fm to a 
value for / of 0.26, which is yet another indication that this regime for the coupling constant 
corresponds to volumes where the instanton effects are important. To obtain this number 
we correct for the finite renormalizations appearing in eq. ( |119| ) and we assumed that setting 



the linear size L of the e- dimensional torus to R gives the MS scheme. This ambiguity could 
be fixed by computing the effective hamiltonian using e.g. Pauli-Villars regularization and 
the known relation to the minimal subtraction scheme. In our calculations / is however 
just a parameter that allows us to probe different volumes and its precise relation to the 
physical scale is not so relevant. 

Returning to fig. |T^, we can distinguish three regimes in R. For small R, we expect 
the finite-size effects, like the effect of the curvature of S 3 , to be large. The masses in the 
effective model, although perturbatively calculable, will therefore not correspond to the 
masses on the torus. In this regime the smallest mass in the 000-odd sector 0~ (see fig. |9[) 
is actually smaller than the mass of the scalar glueball + , an effect that also can be seen 
from the perturbative evaluation of the energy levels eq. ( |163|) and (|164|) . This we expect 
to be a consequence of the strong finite size effects that are visible at small values of R. 

Around R = 0.1 fm, corresponding to L — 0.36 fm, we see (in the sphere geometry) 
the masses deviate from the perturbative expansion, signalling the onset of the instantons. 
They drive the tensor to scalar ratio in the right direction, that is, towards the large- volume 
value of 1.5. The mass of the odd glueball is now higher than that of the even glueball. 
Apparently the non-perturbative effects also drove this quantity towards more physical 
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values: from lattice calculations a value of m(0~)/m(0 + ) ~ 1.7 can be extracted |36 
compatible with what is found here. 

For R > 0.2 fm, our effective model is no longer valid, due to the spreading out of the 
wave functional over the boundary of the fundamental domain, as discussed before. 

Finally, an important goal of this project was to get results for glueball masses as a 
function of 9. The truncated results showed a pronounced dependence on 9 in the regime 
where boundary effects are appreciable (see fig. |6|). Quite remarkably, and unexpectedly, 
this strong dependence disappears when adding the one-loop corrections. In particular, at 
9 = and 9 = it, where boundary conditions can be implemented most accurately, the 
masses do not differ significantly (see fig. |10D . Caution needs to be applied in concluding 
that the same will hold at large volumes, but in any case it would be interesting if glueball 
masses could be measured at 9 = it on the lattice as comparison. 

In conclusion, we should expect our one-loop corrected result to be a relatively accurate 
representation of the true masses on S 3 up to / = 0.28 corresponding to a circumference of 
approximately 1.3 fm, up to where also the variational basis does not exhibit too much of 
the problems with the Temple bound (cf. fig. ^). The approach to infinite- volume values 
of results on the three-sphere, as compared to torus results, is slow. We typically have 
a dependence on powers of 1/R as compared to an exponential behaviour in L ||16|| . As 
our results should not be expected to be already in the asymptotic domain, the rough 
agreement we find with results on the torus is gratifying. 

When comparing the truncated and one-loop corrected results, the results for the 9 
dependence show that strong non-linear and non-perturbative effects influence the spec- 
trum. We have shown there is a small, but finite window from R = 0.1 fm to R = 0.2 fm 
(at smaller volumes everything can be described perturbatively) where these effects can 
be included reliably, showing convincingly how important the global properties of the field 
space are for the non-perturbative dynamics of non-abelian gauge theories. 
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Appendices 

A Details on the basis 

In this appendix we will exploit the relation between S 3 and SU(2) to rederive the ba- 
sis flT2"|). After this we will make the link with the definition of functions on S n in terms of 
homogeneous polynomials in the coordinates of the embedding space. 
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A.l The SU(2) group structure of S 3 

There is a one-to-one correspondence between points of S 3 and SU(2) given by g = n ■ a. 
For instance, the matrix V- of eq. (Q) can in this light be seen as the adjoint representation 
of SU(2). The group of spatial symmetries of S 3 is SO (4). These symmetries correspond 
to the left and right symmetries on SU(2). Let ip be a function on SU(2) or equivalently 
on S 3 . SU(2) can now act on this function by left or right multiplication: 

(v L ( gi )4>)(g)=ij(g^g), (V R (g 2 )^) (g) = ^(gg 2 ), g, 9l , g 2 e SXJ(2). (169) 

These infinite-dimensional representations of the group SU(2) induce representations of the 
Lie algebra su(2). We will show that the generators £ of these representations correspond 
to L\ and L 2 : 

(iCftf))(g) = ^V(exp(-#|r j )^)| t= o 

d i 
= -^ip{n ■ a + ^tajn ■ a)\ t=0 

= ^(( n « + :Nj>^K) lt=o 

= lrf m n,^(n) = (iL{^)(n). (170) 

According to the Peter- Weyl theorem |^7|, a basis of functions on SU(2) is provided by the 
collection of matrix elements of all the finite-dimensional unitary irreducible representations 
of SU(2). The relation of these matrix elements with the basis of eq. (|T2| ) is 

(g\lm L m R ) = (-1)™ L v / 2TTT D l (g) . (171) 



A. 2 The eigenfunctions of the laplacian on S" 

In this section we summarize some well-known facts about the spectrum of the laplacian on 
an n-dimensional sphere. This allows us to identify the functions (O) with explicit func- 
tions on S 3 . Let xo, . . . ,x n be coordinates in IR n+1 . Using radial coordinates r, 9i, . . . , 9 n 
one can explicitly solve for the eigenfunctions of the spherical laplacian A n , using separa- 



tion of variables and recursion in n. The result of this calculation (cf. e.g. [57|) is that the 
spectrum is given by 

{-L{L + n- l)} L=0 ,i,2,.., (172) 
and the degeneracy of the level L is given by 

(173) 

We will relate these eigenfunctions to polynomials in x^. Let V L be the set of polynomials 
in that are homogeneous of degree L. To such a polynomial p corresponds a function 
7 on S" defined by 

p(x , ...,x n )=r L Y(x , ...,x n ). (174) 
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Suppose that p is a harmonic polynomial, i.e. Ap = 0. We then have 



= f^TT^TT + \ A n) (r L Y)=L(L + n- 1) r L ~ 2 Y + r L ~ 2 A n Y. 



(175) 



We have found that every harmonic polynomial of degree L corresponds to a function on the 
n-sphere that is an eigenfunction of the spherical laplacian with eigenvalue —L(L + n — 1). 
The number of monomials of degree L in the variables xq, . . . ,x n is given by 




dim(V L ) = \ J . ( 176) 

Hence we have 

dim ( Ker (A)) = dim (V L ) - dim ( Im (A)) 
> dim (V L ) - dim (V L ~ 2 ) 

L + n\ ( L -2+n 
n I \ n 



(177) 



Since we know the degeneracy of the eigenvalue —L{L + n — l), we have shown that all the 
eigenf unctions of the spherical laplacian correspond to harmonic polynomials. Moreover, 
the inequality in the last equation is actually an equality. 

For the case of S 2 the eigenfunctions are the well-known spherical harmonics Y^. For 
the case S 3 the spectrum becomes {— L{L + 2)}x, = o,i,2,... , with the degeneracy given by 

= {L + lf. (178) 

When comparing this with the results (|171|) , we see that L = 21. The four I = | modes 
(x\ \ mi ran) are linear combinations of the four scalar functions x M , (/x = 0, . . . , 3), whereas 
the nine I = 1 modes {x\1w,l tur) correspond to the nine components of V-. 
In section § the results of this appendix are applied to S 8 . 





B The linear term in the path integral 

In this appendix we show that the term J V Q V in eq. ( |116| ) does not make a contribution 
to the effective action up to the order we are interested in: the lowest order contribution 
from J to the effective potential will be of the form c 4 d 2 and c 2 d 4 . 

The crucial remark is that in the path integral the constant modes are excluded from 
the integration over Qq, and the (c,d) modes are excluded from the integration over Qi. 
This means that one must remember to put projectors (1 — P) around operators like W^, 
and replace source currents like J M with ((1 — P) J) M . We will show that after this projection 
J contains only terms that are cubic in B. 

We have F 0i = -F i0 = B>i and F {j = -2e ijk (c k - d ■ V k ) + [B h Bj], with B i = c i + d-V i . 
This implies Jo = D^F^q = —[Bi, Bi\ and Jj = D^F^j = Bj + D^F^. In our approximation 
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we are only interested in the term -Dj-F^- = J- + J^ 1 + jf \ where is of order n in 
the field B. From = diF^ — e^F^ + \A±, F^], we obtain 

jf ] = -4c i5 (179) 

Jj ^ = 3 Eijk[Ci, Cfc] 3 £rnpn\d m i d^Vj ', (180) 

jj 3) = [Si, [/-?i, .B,-]] 

= [q, [ci, Cj}} + [d ■ Vi, [d ■ Vi, Cj}} 

+[c i ,[c i ,d-V j ]] + [d-V i ,[d-V i ,d-Vj]] 
+ [d ■ Vi, [ci, cj]] + [q, [d ■ Vi, Cj]] 

+[d ■ V, [a, d ■ Vj\] + [a, [d -v t ,d- Vj\\ (181) 

This directly implies that (1 — Py)J^ = (1 — Pv)J ( ' 2 ^ = 0. The first four terms of do 
not contribute since they too are modes with (/, k) = (0, 1) or (/, k) = (1, 0). We write the 
last four terms of as follows: 

A\jVf + B^V^VJ 1 , (182) 

with the su(2)-valued constants 

= [d p , [d, cj]] + [ci, [d p , Cj]], B'1 m = [d m , [c h d n ]] + [a, [d m , d n ]]. (183) 

The matrices A p and Bi can be decomposed in a trace part, an antisymmetric part and a 
symmetric traceless part. For the A p terms, this is the decomposition into (l,k) is (1,0), 
(1,1) and (1,2) respectively. For the Bi terms, it gives the split in the three possible I 
values (0, 1), (1, 1) and (2, 1). Projecting out the (c, d) modes from Jj, we obtain 

((1 - P v )J)j = A%VP + BrVTVp, (184) 



where A and B are obtained from the A and B of eq. ( |183j ) by projecting out the trace 
part. 

Expanding the path integral of eq. (|116[) gives diagrams with interactions of the Q field 



with the source J{B). If we also take along the Q 3 , BQ 3 and Q A terms in the action, 
diagrams with only one occurrence of J can in principle give a contribution of the order 
c 2 d 2 . However, in order to obtain a c 2 d 2 invariant, we must take the trace part of A and 
B. Since these trace parts are zero, we see that a single occurence of J can not give rise 
to contributions of the order c 2 d 2 . 



C Angular matrix elements 



We will perform the reduction of the matrix elements of the operator A of eq. ( |145|) to 
reduced matrix elements. We use eq. ( |132| ). 



(j , ,m , ,i' 1 ,i , 2 \A\j,m,i 1 ,i 2 ) 

{|0, 0, l s ; Li, 0, fi; L 2 , 0, f 2 ) \j, m, l s ; L x , l x , n; L 2 , l 2 , r 2 )} 
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E (-i) 



2' 1 ,m' 2 



E (-i)^ m: 



mi ,m2 



2/ + 1 

27+1 
i 



i'i i' 2 f 



x 



i "<>2 
h h j 

m l m 2 — m 



m 



x 



m',m s ,ni s 



-1) 



1 



-1) 



S " b S 



I2V S + 1 ' sj2l s + 1 ' V2/ s + 1 

(I^; Z' s , r(; m' s , m'^Li] l s , 0, n; m s , 0)|Li; i s , /i, n; m s , mi)} x 
(4; C Z 2> r 2! - m ' s > m ' 2 l{l^2; 0, f 2 ; -m s , 0)|L 2 ; i„ Z 2 , r 2 ; -m s , m 2 )}. 

We introduce reduced matrix elements through 

(L[; l' s , l[, t[; m' s , mi|{|Li; l s , 0, fi; m s , 0) |Li; Z„ /i, n; m s , mi)} = 



(185) 



■1) 



i' s +m' 1 



I' 



-m' m, m. 



l[ h 
-m' x mi 



F(ii,ii,zi). 



(186) 



Using 



Z'i h 
-m' l mi 



■1) 



ll—mi r r 



(187) 



the summations over m! x and m' 2 become trivial. The remaining summations in eq. (|185|) 
can then be dealt with using twice the completeness relation 



£ | 3i 32 33 

mi m 2 \ m l m 2 m 3 



3i 32 3 3 
rrii m 2 m 3 



2j 3 



S (3lJ2,33), 



(188) 



where S(ji, j 2 , j 3 ) = 1 if ji, j 2 and j 3 satisfy the triangular condition, and is zero otherwise. 
We obtain 



(j',m f , i' 1 ,i' 2 \A\j,m,ii,i2) = 



F(i'i,ii,ii)F(i' 2 ,i 2 ,i 2 ) 
(2/^ + l)(2/; + l)(2/ s + l) 



(189) 



We absorbed the factors (—1)^(2^ + 1) for i — 1, 2 in the F functions. We can evaluate 
the F functions by setting in eq. (|186|) 



mi = h, m s = l s , m' = I' 



m. 



I' - Is- 



(190) 



V. 



Is 



Is 



-1 



Thus we find 

F(i',i,i) = (-l) l '° 

(L'\ l' s , l r , r ; l' s , l' r \ \ \L; l s , 0, f; l' s — l s , 0) \L; l S7 l r , r; l s , l r 



-I' V - I I 



(191) 
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Note that the only dependence of the matrix elements on j is through the triangular 
condition on Zi, l 2 and j. Also note that the explicit S functions 

6 lih 6 Vah 6(l'J a ,l a ) (192) 

are superfluous: they are also contained within the F functions. The 5(li,l 2 ,j) function 
can be deleted too: the corresponding triangular condition is satisfied from the beginning. 

Specializing to the case where A only depends on c, we have l s = and we can also 
absorb the factors (—l) la (2l s + 1)~ 1//2 in F. We have for instance 

(j , ,m',i' 1 ,i , 2 \mv 3 (c)\j,m,i 1 ,i 2 ) = 5 j > j 5 m , m F 3 (i' 1 ,i 1 )5(i 2 ,i 2 ), (193) 

with 

F 3 (i',i) = (L';l' s ,l' r ,T';l' s ,l' r \mv 3 (c)\L;l s ,l r ,T;l s ,l r ). (194) 

We similarly treat inv 4 (c), inv 6 (c), inv^c) and inv 8 (c) which give rise to the reduced matrix 
elements F 4 , F e , F 7 and F s respectively. Operators that depend only on d and operators 
that are products of a function of c and of a function of d pose no problems either. We 
have for instance 

(j',m',i' 1 ,i' 2 \iav 3 (d)\j,m,i 1 ,i 2 ) = 5 j i j 8 m r m 8(i' 1 ,i 1 )F 3 (i 2 ,i 2 ), (195) 
(j\ m\i[,i' 2 \ my 3 (c)iav 3 (d)\j, m,!!,^) = 8j>j8 m < m F 3 {i' l ,i x )F 3 {i 2 ,i 2 ). (196) 

The operator mv 2 (c,d) is treated with the general formulae eq. ( |189| ) and ( |191[ ). One can 
show inv 2 (c, d) = § - ±V5(cd\0, 0, 2; 2, 0, 1; 2, 0, 1). Using eq. (PI we define 

iY* 

(F; l' s , l' r , t>; l' s , Q {|2; 2, 0, 1; l' s - l s , 0)\L; l s , l r , r; l s , l r )} , (197) 
and we obtain 

0", m', i[, i 2 \ inv 2 (c, d) | j, m, i u i 2 ) = 

1 ^22(^^1)^22(^2^2) 

[(2j; + i)(2/. + i)p 

For products like inv 3 (c) inv 2 (c, d) we obtain 

0",m', i'^i' 2 \ inv 3 (c) mv 2 (c, d)|j, m, ii, i 2 ) = 

r r J 1 771 /- •/ • \ r/ •/ • \ 1 -^52 (jj , ^l)-p22 (^2 jjg) j /-, 00 s 

OfjO m > m < gF 3 [i u H)o{i 2 , % 2 ) - 997777 77777 77777 \ , (199) 



x r I lr/-/ ■ U /.i • \ 1 ^22^1^1)^22(^2^2) I / ln0 N 
Oj'jflm'm i g0{t 1 ,ll)d(l 2 ,l 2 ) ~ 997777 7777 777T } ■ (198) 



[(2/; + l)(2Z a + l)]3 



with 

-1 



6 C fc c fcc to 



(L'; l' s , l' r , t'- C /:| {inv 3 (c)|2; 2, 0, 1; l' s - l s , 0)|L; Z s , Z r , r; Z s , Z r >} . (200) 
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Similarly the operator inv 4 (c) inv2(c, d) leads to the reduced matrix element F 62 . 

The only operator in V 2 that we have yet to deal with is invite, d). For this operator 
we need to construct the spectral decomposition. We write 

inv 2 , = a|0,0,4;4,0,l;4,0,l) + 6|0,0,2;4,0,l;4,0,l) 
+c |0, 0, 0; 4, 0, 1; 4, 0, 1) + d |0, 0, 2; 2, 0, 1; 2, 0, 1) 
+e (|0, 0, 2; 4, 0, 1; 2, 0, 1) + |0, 0, 2; 2, 0, 1; 4, 0, 1)) 
+/ (|0, 0, 0; 4, 0, 1; 0, 0, 1) + |0, 0, 0; 0, 0, 1; 4, 0, 1)) 

+g |0, 0,0; 0,0,1; 0,0,1). (201) 

Using the explicit formulae for the corresponding polynomials, we solve for the coefficients 
a, . . . ,g. The new reduced matrix elements we need are F44, F 42 and _F 40 for the operators 
with (L, l s , l r ) respectively (4, 4, 0), (4, 2, 0) and (4, 0, 0). The function F 40 can be obtained 
from F4 using the relation 

inv 4 (e) = ± + h{i ^ °> °> °> °>- ( 202 ) 

The result is 

inv^c, d)|j,m,ii,i 2 ) = 8j>j8 m > m \ r + B \ , (203) 



with 



7-F 4 4(i / 1 ,ii)F 44 (i 2 ,i 2 ) 



11583 
4 

'50193" 42 ^l' " l >* l - v '-'- '-' 1.859 
1 



^^0^^42(^1,^1)^42(4^2) - 7^0^22(^1,^1)^22(4^2) 



1859 V 2187"( F42 ( i i' ^) F 22(4, ^2) + F 22 (i' 1 ,i 1 )F 42 (if 2 , i 2 )), (204) 
B = ^F 4 (i[, ii)F 4 (i 2 , i 2 ) - ^(F 4 (i[, i 1 )5(i' 2 , i 2 ) + 5(i[, ii)F 4 (i' 2 , i 2 )) 

+^5(4,^)5(4,^). (205) 
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-3-2-10 1 2 

u 

Figure 1: The (u,v) plane: a two-dimensional subspace of A. Location of the classical 
vacua, sphalerons, lines of equal potential. The boundary of the fundamental domain 
lies between the Gribov horizon (fat sections) and the lower bound A (drawn parabola): 
A C A C fl, with Q the Gribov region. 
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Figure 2: Expansion of V^\u) in the tunnelling parameter u. We dropped the e and log(L) 
dependent parts. We have drawn the expansion up to order u 4 , u 6 , u s and u 10 . Longer 
dashes correspond to higher order expansions. The horizontal line at u — 1 denotes the 
exact potential at the sphaleron, the drawn curve is the sixth order fit explained in the 
text. 
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Figure 3: Lowest energy levels for 6 = 0. Drawn curves correspond to levels in the (0, 0, 0) 
sector. The dashed curve denotes the ground level in the (1, l,2)-even sector. The short- 
dashed curves are the perturbative expansions, and the individual dots are lower bounds 
on the levels as obtained by Temple's inequality. 
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Figure 4: Glueball masses for 6 = as a function of the coupling constant. The drawn 
curves are the masses of the first scalar (0 + ) and tensor (2 + ) glueball. The dashed lines 
denote the lower bounds, the dotted lines the perturbative results. 
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Figure 5: Strong coupling limit, 6 = 0: the drawn and dashed lines are the variational 
results for the lowest levels in the sectors 000-even and 112-even respectively. The dotted 
lines represent the analytic strong coupling limit. 
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m(0 + ) R 
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Figure 6: Scalar glueball mass at / = 0.25, / = 0.3, / = 0.4 and / = 0.5 as a function of 
6. The dashes at 9 — ±7r denote the variational result with exact implementation of the 
boundary conditions. Note that the vertical range in these four plots is the same. 
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Figure 7: The wave function \1> oc \ip(r c , rd)\/(r^rj) plotted at r d = \/3 for the ground state 
in the sector 000-even. 
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Figure 8: One-loop results. Lowest energy levels for 9 = 0. Drawn curves correspond to 
levels in the (0, 0, 0) sector. The dashed curve denotes the ground level in the (1, 1, 2)-even 
sector, the short-dashed curves are the perturbative expansions. The four small plots also 
show the lower bounds (individual dots) as obtained by Temple's inequality. 
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Figure 9: One-loop results. Glueball masses for 6 = as a function of the coupling 
constant. The drawn curves are the masses of the various scalar and tensor glueballs. The 
dotted lines denote the perturbative results. 
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Figure 10: One-loop results. Scalar glueball mass at / = 0.15, / = 0.2, / = 0.25 and 
/ = 0.3 as a function of 9. The dashes at 9 = ±7r denote the variational result with exact 
implementation of the boundary conditions. The vertical range of the plots is the same as 
in fig. ||. 
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Figure 11: One-loop results. The wave function \I> oc \ip(r c , r d )\/(rjr^) plotted at r d = \/3 
for the ground state in the sector 000-even. 
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Figure 12: Comparison of the tensor to scalar mass ratio m 2 +/m + as obtained on the 
three-sphere to the relevant ratios on the torus. The dotted line denotes the perturbative 
expansion, the bar on the right indicates the range of lattice Monte Carlo values at L — 1 
to L = 1.5 fm. 
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